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Abstract We review some recent developments in charm meson physics. In particular 
discuss theoretical predictions and experimental measurements of charmed meson decays to lep- 
tonic, semileptonic, and hadronic final states and implications of such measurements to searches 
for new physics. We discuss D° — D°-mixing and CP-violation in charm, and discuss future 



experimental prospects and theoretical challenges in this area. 
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1 INTRODUCTION 



The charm quark has played a unique role in particle physics for more than three 
decades. Its discovery by itself was an important validation of the Standard 
Model (SM), as its existence and mass scale were predicted ([I]) on the basis of 
low-energy kaon experiments before any direct experimental signature for charm 
was available. 

Several features distinguish charmed hadrons from those with other flavors. 
While their mass, 0(2 GeV) places them in the region where non-perturbative 
hadronic physics is operative, theoretical methods developed for heavy quarks 
can in principle still be applied, albeit with larger uncertainties. On the other 
hand, recent advances in unquenched lattice QCD simulations paved the way for 
charm data to be used to probe the Yukawa sector of the SM. It is the only 
up-type quark that can have flavor oscillations. Finally, charm decays provide a 
unique window on new physics (NP) affecting the up-type quark dynamics. In 
many cases, charm transitions provide almost background-free low-energy signals 
of NP. For example, many popular NP models predict signals for CP-violation 
(CPV) much larger than what is generally predicted within the SM (PJ). One 
hopes that, just like manifestation for charm quark existence came from low- 
energy kaon oscillation experiments, oscillations of charmed hadrons can provide 
hints of what is happening at the TeV scale. 

Experiments operating at the ^(3770) resonance, near threshold for DD pro- 
duction, such as MARK III at SPEAR, performed the initial exploration of charm 
phenomenology ([3]). Later, higher energy machines, either fixed target exper- 
iments operating at hadron machines or higher energy e + e~ colliders, entered 
this arena, with much bigger data samples. In recent years, we have seen a re- 
newed interest in studying open charm in e + e _ colliders with a center-of-mass 
energy close to DD threshold. The CLEO experiment (.4) at CESR, after years 
of charm studies at a center-of-mass energy near the T(45), has collected a sam- 
ple exceeding 800 pb^ 1 at the -0(3770) center-of-mass energy, and a sample of 
about 600 pb _1 at a center-of-mass energy close to 4170 MeV, optimal for D s 
studies. The BES-II experiment, at BEPC, has published results based on 33 
pb _1 accumulated around the ^(3770). It has an upgrade program both for the 
detector (BESIII) and the machine (BEPCII), designed as a charm factory with 
10 33 cm _2 s _1 peak luminosity that should be completed in 2008 ([5]). In parallel, 
BABAR and Belle have provided important contributions to our knowledge of 
charm decays, exploiting their impressive data sets. Both B factory facilities 
have achieved luminosities almost ten times their original design. A KEK-B up- 
grade has been approved to attain luminosities about one order of magnitude 
greater than already achieved. An alternative approach, with higher luminosity 
goals and the added capability to run at both T(AS) and at the ^(3770), is being 
considered, but is not approved. Finally, CDF and DO have entered the arena 
of charm physics, applying to this study some of the tools developed to pursue 
beauty physics. 

Experiments operating at the ip(3770) resonance have several advantages: the 
final state is extremely simple, being dominated by a DD pair. The cross section 
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for charm production is relatively high, a(D°D°) = 3.66 ± 0.03 ± 0.06 nb and 
a(D + D~) = 2.91 ± 0.03 ± 0.05 nb ((6j). In addition, the relatively high branching 
fractions to low multiplicity final states allow the use of tagged samples, where 
one D is fully reconstructed (tag) and the rest of the event recoiling against 
the tag is used to study specific signals. This technique is particularly useful to 
study leptonic and semileptonic events, as it allows a precise reconstruction of the 
neutrino 4-momentum, and of the invariant mass squared of the lepton-neutrino 
pair (q 2 )- In addition, the DD pairs are produced in a C = — 1 state, and this 
quantum coherence allows unique probes of mixing and CPV . On the other 
hand the b-factory experiments, Belle and BaBar, can also be considered charm 
factories. At the T(4S) center-of-mass energy, a(bb) ~ 1.1 nb, while ct(cc) ~ 1.3 
nb. The background to be tackled is higher than at the lower energy, and neutrino 
and q 2 reconstruction in general are not as precise. Significant improvements are 
achieved, however, through the use of tagged samples ([7]), made possible by the 
impressive size of the data set accumulated: in January 2008 the samples are 484 
fb _1 at BABAR and 763 fb _1 at Belle. Also, the energy at which they operate 
makes possible the production of charmed baryons and boosts the charm particles 
sufficiently so that time dependent measurements are possible. 

Experiments at hadron machines have the advantage of much higher cross sec- 
tions for charm and beauty production, at the expense of significant backgrounds. 
Luckily the relative long lifetime of charm hadrons (~ 1 ps), combined with the 
development of silicon micropattern detectors provides a unique and powerful 
signature of charm meson decays: the identification of detached secondary ver- 
tices. About 30 years after its initial inception ([8]), this technique is still being 
perfected, introducing vertex algorithms more and more sophisticated that pro- 
vide charm and beauty event tagging almost in real time. This development has 
allowed fixed target experiments and the two Tevatron experiments, CDF and 
DO, to provide significant contributions to our knowledge of charm decays, with 
competitive limits on some rare decays and recent results in D°D mixing. This 
work will be continued by LHCb, the first dedicated charm and beauty experi- 
ment at a hadron collider (LHC), that relies heavily on detached vertex criteria 
in early stages of its triggering process and is considering an upgrade that would 
include detached vertex criteria in the lowest trigger level ([9]). 



2 LEPTONIC AND SEMILEPTONIC DECAYS 

Charm leptonic and semileptonic decays are ideal laboratories to study non- 
perturbative QCD, and to determine important quark mixing parameters. In 
addition, they may provide additional constraints on physics beyond the SM. 
In the SM semileptonic decays are described by an effective Hamiltonian 

Gf 

H = -^V cq L»qT^ (1) 

with being a leptonic current, Gf is a Fermi constant, and T M = 7^ (1 — 75). 
Theoretically, leptonic decays are the simplest to describe, as they only depend 
on a single non-perturbative parameter, the decay constant fo q , 

(0\q^ l5 c\D q ) = -if D y D , (2) 

which parameterizes the amplitude of probability for the heavy and a light quark 
to "find each other" in a meson. Semileptonic decays are traditionally described 
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in terms of form factors cast as a function of q 2 , the invariant mass of the electron- 
neutrino pair. Experimental determinations of these form factors are performed 
through the study of the differential decay width dT/dq 2 . In both cases, decay 
constants and form factors are QCD parameters that can only be computed using 
non-perturbative techniques. 

Accurate calculations of non-perturbative QCD parameters are very challeng- 
ing. Lattice QCD represents an appealing approach. In principle, it is the only 
one that can be improved in a systematic way. A big stumbling block has been the 
inclusion of dynamical quark effects (unquenched lattice QCD). Recently, tech- 
nical developments such as highly improved actions of QCD and the availability 
of "2+1 flavor" MILC configurations with 3 flavors of improved staggered quarks 
have lead to results with much higher accuracy and allowed for consistent esti- 
mate of both statistical and systematic errors involved in the simulations. Two 
groups have reported charm decay constant calculations with three dynamical 
quark flavors: the Fermilab/MILC Lattice collaboration (10) and the HPQCD 
collaboration (jTTj) . They both use the "2+1" MILC configurations including three 
flavors of staggered quarks: one heavier with a mass close to the strange quark 
mass m s , and two degenerated light quarks with masses between m s /10 and m s , 
but differ on how they treat heavy quarks in their formulations of lattice QCD. 
The Fermilab group has also calculated the shape and normalization of the form 
factors in semileptonic D — > Trev e and D — » Kev e decays f|12|) . 

The QCD sum rules fll3|.ll4|) provide a method for calculating hadronic matrix 
elements including non-perturbative effects that was designed to make maximum 
use of the known manifestations of non-perturbative QCD. A few parameters 
describe the non-pertubative dynamics and are fixed from well known hadronic 
processes, and then applied to heavy meson decays. Finally, quark models, gen- 
erally QCD inspired and based on a variety of assumptions, have been used to 
predict form factor normalizations and decay constants (|16H . 

In principle, charm meson semileptonic decays provide the simplest way to 
determine the magnitude of quark mixing parameters: the charm sector allows 
direct access to \V CS \ and \V c d\- Semileptonic decay rates are related to \V cq \ 2 via 
matrix elements that describe strong interaction effects. 

The study of charm semileptonic decays may contribute to a precise deter- 
mination of the Cabibbo Kobayashi Maskawa {CKM) matrix element \V U \,\. A 
variety of theoretical approaches have been proposed to use constraints provided 
by charm decays to reduce the model dependence in the extraction of \V u b\ from 
exclusive charmless B semileptonic decays. In particular, if heavy quark effective 
theory (HQET) (fTT|) is applicable both to the c and b quarks, there is an SU(2) 
flavor symmetry that relates the form factors in D and B semileptonic decays 
(|18j) . For example, a flavor symmetry relates the form factors in D — > itlv are 
related to the ones in B — > ttIu, at the same E = v ■ p n , where E is the energy 
of the light meson in the center-of-mass D frame, v is the four-velocity of the D 
meson, and pp is the 4- momentum of the light hadron. The original method has 
been further refined (|19j) : the large statistics needed to implement these methods 
may be available in the near future. 
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2.1 Theoretical Predictions for the Decay Constant 

The leptonic decay width is given by 

- iv) = %fh q m\M Dq (l - ^pj \V cq \ 2 , (3) 

where q = d,s for D + or D s states respectively, Mu is the D q mass, mi is the 
mass of the final state lepton, and \ V cq \ is the CKM matrix element associated 
with the c — ► q transition. Due to helicity suppression, the rate goes as mf; 
consequently the electron mode D + — > e + v e has a very small rate in the SM. The 
relative widths scale as 2.65 : 1 : 2.3 x 10 -5 for the t + z/ t , /i + z/^ and e + u e final 
states, respectively. The decay constant fo parameterizes a matrix element of 
the axial current and is the only non-perturbative parameter in Eq. (|3|). It can be 
related to the wave function overlap of charm quark and light antiquark. Charm 
meson decay constants are amenable to experimental and theory determination 
both for D + and D s , thus allowing a direct measurement of SU(3) breaking and 
a comparison with the theory. 



2.2 Experimental Determinations of fjj 

The CLEO collaboration (27) has measured f D + = (222.6 ± 16.7±|;|) MeV, using 
a tagged sample of D + D~ decays collected at a center-of-mass energy close to 
3.77 GeV. The existence of the neutrino is inferred by requiring the missing mass 
squared (MM 2 ) to be consistent with zero, 

MM 2 = (E beam - E^ + ) 2 - {p D - -p^) 2 , 

Figure Q] shows the measured MM 2 , with a 50 event peak in the interval [-0.050 
GeV 2 , +0.050 GeV 2 ], approximately ±2a wide. The background is evaluated as 
2.81 ± 0.30 db 0.27 events. The same tag sample is used to search for D + — > e + u e . 
No signal is found, corresponding to a 90% CL upper limit B(D + — > e + z/ e ) < 
2.4 x 10 -5 . More data is available on f D +- Early measurements and a recent 
BaBar result of f D + determine the ratio B(D S — * fj,v)/B(D s — > (f>ir) (|28p . This 
adds an additional large source of error as the denominator is not well known (|33p . 
CLEO-c uses a sample of D S D* collected near the center-of-mass energy of 4.17 
GeV to study D s leptonic decays (34). They measure the branching fraction for 
the decays D s — > /iv^, D s — > tu t , with r — > irv T , and r — ► ev T v e ([35]) . Recently, 
Belle has reported an absolute value for the branching fraction B(D S — > (J,Vu) 
based on fully reconstructed samples of events of the type e + e~ — > D*DKX, 
D* — > Dsj where X is any number of tt and at most one 7 from fragmentation 
(|36p . Theoretical predictions are summarized in Table [lj while measurements are 
shown in Table [2j The average of the absolute measurements is f n + = 275 ± 10 
MeV, assuming that |V^ S | = \y c d\ = 0.9737. Typically, the experimental value is 
above theoretical predictions. In general, the errors are such that the discrepancy 
is not yet meaningful, with the exception of the most recent UKQCD-MILC 
calculation pip . In this case, the discrepancy between theory and experiment 
exceeds their stated errors by about 3 a. 
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2.3 Constraints on New Physics from f D 

Leptonic decays are sensitive probes of NP interactions mediated by charged 
particles. Models with an extended Higgs sector, which include new charged 
scalar states, or models with broken left-right symmetry, which include massive 
vector Wft states, are primary examples of such interactions. Recent evidence 
of observation of B — > tv t decay brought renewed attention to such models. In 
particular, two Higgs doublet models, including Minimal Supersymmetric SM 
(MSSM), could give contributions to such transitions. Different implementations 
of this extension of the SM can be formulated (I29p . For example, the first doublet 
(3>i) could give mass to the up-type fermions and the second ($2) to the down- 
type fermions. In this case, 

B(D+^£+u e ) = B SM 

B(Df^£ + u e ) = B SM 

Note that the latter model introduces a correction to the SM expectations that 
may be considerable and negative at large tan 2 (3. A limit can also be set on 
the mass of a charged Higgs, m H + > 2.2 tan (3. This limit is similar to the one 
obtained from the measurement B(B — > rv) decay (|30p . 
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2.4 Absolute Branching Fractions for Semileptonic D Decays 

Determination of absolute branching fractions for D semileptonic decays consti- 
tute important measurements. Assuming \ V CX \ to be known, they determine form 
factor normalization. Conversely, if the form factors are known independently, 
say, from the lattice QCD calculations, these branching determine the relevant 
CKM matrix elements. By comparing the inclusive branching fractions of the 
D + and D° mesons with the sum of the measured exclusive branching fractions, 
one can determine whether there are semileptonic decay modes as yet unobserved. 
BES-II( [52l) and CLEO-c (l3ip have recently presented data on exclusive semilep- 
tonic branching fractions. BES-II results are based on 33 pb _1 ; CLEO-c's results 
are based on the first 57 pb _1 data set. Both experiments use tagged samples 
and select a specific final state through the kinematic variable: 

U = E m i ss |cp m j ss |, 

where E m i ss represents the missing energy and p represents the missing momen- 
tum of the D meson decaying semileptonically. For signal events, U is expected 
to be 0, while other semileptonic decays peak in different regions. Fig. [2] shows 
the U distribution for 5 exclusive D + decay modes reported by CLEO-c, which 
demonstrate that U resolution is excellent, thus allowing a full separation be- 
tween Cabibbo suppressed and Cabibbo favored modes. Table [3] summarizes the 
recent measurements from CLEO-c and BES-II, as well world averages reported 
in the Review of Particle Physics (33). 

Absolute branching fractions for D° — > Klu have been recently published by 
Belle ©: they obtain B(D° -> Klu) = (3.45 ± 0.07 ± 0.20)% and B(D° -» 
tt£u) = (0.255 ± 0.019 ± 0.016)%. CLEO-c uses the two tagging modes with 
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lowest background (D — > i"T + 7r and D — > K + ir it ) to measure the inclusive 
D° and D + semileptonic branching fractions (|42h . They obtain 

B{D + le+^e) = (16.13 ± 0.20 stat ± 0.33 sys )%, 

£(D° -> Xe + i/ e ) = (6.46 ± 0.17 stat ± 0.13 sys )%. 

The sum of the exclusive semileptonic absolute branching fraction is B{D + — ► 
Xe + u e ) exci = (15.1 ± 0.5 ± 0.5)% and 5(L>° -> Xe + ^) cxd = (6.1 ± 0.2 ± 0.2)%: 
the measured exclusive modes are consistent with saturating the inclusive widths, 
although there is some room left for higher multiplicity modes. The CLEO-c data 
have been used in this comparison, as they dominate the present world average: 
the exclusive modes are consistent with saturating the inclusive semileptonic 
branching fraction at a 41% confidence level in the case of the D + and 18% 
confidence level in the case of the D°. 



2.5 Form Factors For The Decays D — > K(tt)£v 

Theoretical parameterizations of semileptonic decays involve two non-perturbative 
quantities parameterizing matrix element of a single hadronic current. Tradition- 
ally, the hadronic matrix elements for transitions to pseudoscalar hadrons are 
described in terms of two form factors, f+(q 2 ) and f-(q 2 ), 

(K(n)\qT^c\D) = f + {q 2 )P» + f-(q 2 )q», (5) 

where P = po + Pk(w) an d q = Pd ~Pkm- An alternative parameterization is 
also often used, 

(m 2 n — m\r, \ \ m 2 n — m/L, n 

p*- q 2 ( y ]/+(g 2 )+ f Kw <ffotf), (6) 

with fo{q 2 ) = f+{q 2 ) + f-{q 2 )q 2 /{ m2 D — m x(7r))' Experimental determinations 
of these form factors are performed through the study of the differential decay 
width dT/dq 2 . For cases where the lepton in the final state is an electron and has 
a negligible mass with respect to the parent D, only a single form factor, f + (q 2 ), 
contributes. The partial decay width is given by 

dT(D^K(7r)eu e ) _ G 2 F \V cq \ 2 3 2 2 

df " 24tt3 p ^)\U{q )l (7) 

where PkM the hadron momentum in the D rest frame. Form factors have been 
evaluated at specific q 2 points in a variety of phenomenological models (|43h . where 
the shape is typically assumed from some model arguments. In order to restrict 
the function space studied, a dispersive representation (j44|) allows to place rather 
general constraints on the shapes of the form factors from their analytic prop- 
erties. Particular parameterizations of the form factors are nevertheless useful. 
The most common parametrization has been a single pole form factor, where the 
pole is the lowest mass resonance formed by the initial and final state hadron. 
For example, in the decay D — ► ireve the dominant pole is the D* . Now that 
more precise data are available, more complex representations are investigated. 
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One class of parameterizations include the dominant pole form factor and ap- 
proximates the dispersion integral by a number of effective poles 

(1 — a) 1 — (q /mf- ' ' 1 ' '.' 



where a determines the strength of the dominant pole, pk gives the strength of the 
kth term in the expansion, and 7& = m v /m v , with my k representing masses of 
the higher order poles. The true form factor can be approximated to any desired 
accuracy by introducing a large number of finely spaced effective poles. In effect, 
it is desirable to keep the number of terms in this expansion to a manageable 
number. The popular Becirevic-Kaidalov (BK) parametrization (|45p is a simpli- 
fied version of the N = 1 truncation of this expansion. In general, both the N = 
case (simple pole) and the N = 1 case can provide good representation of the 
data if the pole masses are allowed to be non physical. An alternative approach 
(150p utilizes a series expansion around an arbitrary q 2 value to- To achieve a 
convergent series, the expansion is formulated as an analytic continuation of the 
form factors in the complex t = q 2 plane. There is a branch cut on the real axis 
for t > (Mp + Mk,-k) 2 , which corresponds to a region associated with production 
of states with appropriate quantum numbers. The convergence is accelerated by 
mapping the whole cut region onto the unit disk z < 1, where z is defined as 



z(q ,t ) = -j==— -==, (9) 

where t± = {Mo ± M^ j7r ) 2 and to is the arbitrary g 2 value that maps onto z =0. 
The form factors are then expressed as (|5Up 

^ oo 

Uiq2) = P(q 2 )Hq 2 ,t ) J] ak{t0Mq2,t0)]k (10) 

with P(q 2 ) = z(q 2 ,my), which accounts for the pole in the form factor at q 2 = 
m v . The physical observables are not expected to depend on $(g 2 ,to) ; which 
can be any analytical function, or to- 

Unquenched lattice QCD calculations for D — > Klv and D — > irlv have 
recently been reported (|12p. The chiral extrapolation is performed at fixed 
E = v-pKfay The lattice "data points" are fitted to the BK parametrization (|45p 

f + (q 2 ) = n „ 2V (11) 

(1 - q 2 )(l - aq 2 ) 



fo(q 



2^ /+(0) 



l-q 2 /^ 



where q 2 = q 2 /m 2 n ,, and a and j3 are fit parameters. The fitted parameters 
are shown in Table HI The FOCUS experiment (46) was the first to perform 
a non-parametric measurement of the shape of the form factor in D — ► KjjLU^ 
(pxLD . CLEO-c (|52]), Belle ©, and BaBar fl53|) reported similar analyses. Fig. [3] 
shows the lattice QCD predictions for D — > Ktv and D — ► with the Belle 
data points superimposed. Table [5] summarizes the experimental form factor fits 
compared to the lattice QCD predictions. By combining the information of the 
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measured leptonic and semileptonic width, a ratio independent of \V c d\ can be 
evaluated, which can serve as a check of the theoretical calculations. For instance, 
assuming isospin symmetry, i.e. F(D — ► ire + v e ) = T(D° — ► 7r _ e + ^ e ) = 2r(D + — > 



7r°e + i/ e ), a ratio 



i? = ^T(D+ /ii/ M )/r(I> 7r e +i/ e ) 

can be formed. Using the recent unquenched lattice QCD calculations (|!Ul ll2 p . 
this ratio can be computed to be 



The quoted error is evaluated through a careful study of the theory statistical 
and systematic uncertainties, assuming Gaussian errors. The corresponding "ex- 
perimental" ration can be calculated using the CLEO-c fu and isospin averaged 



K7 = \ I r exp (-D + -> v-v) _ 249 ± 022 

V T ex P(D -> irev e ) 



The theoretical calculations and data are consistent at 28% confidence level. 



2.6 The CKM Matrix 

An important goal of the next generation of precision experiments is to perform 
direct measurements of each individual parameter. This will enable us to perform 
additional unitarity checks with precision similar to that achieved currently with 
the first row (|47p . With the help of the unitarity constraints, charm-quark-related 
V c d and Vcs are now determined with rather high precision (33). The most recent 
results from LEP II, using the W — > Iv branching fraction, and additional inputs 
from other CKM parameter measurement is V cs = 0.976 ± 0.014 (j48j) . The 
unitarity constraint implies V cd ~ V us = 0.2227 ± 0.0017 (|47j) . 

CLEO-c (49) has extracted \V c d\ and |V CS | by combining the |V^ 9 |/+(0) results 
from the three parameter series expansion fit (|50p with the unquenched lattice 
QCD predictions for /+(0) CG]) to obtain 

\Vcs\ = 1.015 ±0.010 ±0.011 ±0.106 
\V cd \ = 0.217 ±0.009 ±0.004 ±0.023 

The first two errors are experimental, statistical and systematic, while the last 
errors are theoretical, dominated by the discretization uncertainties in the lattice 
QCD charm quark action, which should be improved in the near future. It will 
be interesting to see a unitarity check performed on the second row of he CKM 
matrix element once these errors are further reduced. 



2.7 Form Factors in Semileptonic D — > Vlv Decays 

The structure of the hadronic current in semileptonic decays including vector 
mesons in the final state is more complex, involving four independent form factors, 
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V, Aq, A u and A 2 



(K*(p)\q^c\D) 



2 V ^ 




m D + m K * {p) 



(K*(p)\q W c\D) 



= i {m D + m K , {p) ) e* 



m D + m K *( p ) 




M A 2 {q 2 ) 



(12) 



+ 2im D — —q^ A (q 2 ), 



where e* is a polarization of the final state meson. The vector form factor V 
is dominated by vector meson resonance exchanges, A® is dominated by pseu- 
doscalar meson resonance exchanges, A\ and A 2 are dominated by axial meson 
resonance exchanges. Generally a single pole form factor is assumed, both in ex- 
perimental studies and in theoretical calculations of the normalization of the form 
factors. The FOCUS experiment has developed an interesting technique that ex- 
tends its non-parametric determination of the form factors in D° —* K~/j, + Vn 
to D + — » K~ir + e + v e (55), later adopted also by the CLEO-c experiment (|56|) . 
This method allows model independent determinations of the form factors and 
thus provides a better check for theoretical calculations. For example, Fajfer and 
Kamenic have studied these decays by including contributions of charm meson 
resonances beyond the simple pole (|57p and have found that including two poles 
in the vector form factor improves the agreement between their predictions and 
experimental results. 

The FOCUS experiment (58) has reported evidence for the presence of a small 
even K~ir + amplitude interfering with the dominant K*° component in the decay 
D + — ► K~ir + CLEO-c ([56]) has seen the same effect in D + — > K^7r + e + u e . 
This observation opens up new areas of investigation in exclusive charm semilep- 
tonic decays, namely the investigation of light quark spectroscopy. For example, 
it would be interesting to verify whether this broad s-wave resonance can be iden- 
tified with the k seen in D + Dalitz plot analyses (59). It will be interesting to 
search for similar interference phenomena in D s semileptonic decays. 

3 RARE AND RADIATIVE DECAYS 
3.1 Theoretical Motivation 

Rare charm decays hold great potential to be a sensitive probe of NP. Among all 
rare charm transitions, the most interesting are the decays that are associated 
with AC = 1 flavor-changing neutral currents (FCNC), i.e. transitions that 
change charm quark quantum number by one unit while conserving the electrical 
charge of participating quarks. Examples of such transitions include (a) rare 
radiative decays mediated by c — > wy or c — > wyj quark currents, or (b) rare 
leptonic and semileptonic decays mediated by c — > uii quark currents. Here t 
could either be a charged lepton such as e or [i or a neutrino v. In addition, fully 
non-leptonic FCNC, such as c — ► ug or c — ► uqq are possible. We shall discuss 
them in Section SJ 

In the SM, where FCNC cannot occur at the tree level, this is usually associated 
with large contribution of top quark to one-loop electroweak diagrams due to 
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the Glashow-Iliopoulos-Maiani (GIM) mechanism (|60p . This assures that the 
bottom-type FCNC decay is dominated by the short-distance contributions and 
therefore is reliably computable. It has become evident that this situation is 
not realized in charm decays due to relatively small mass of the bottom quark 
and significant hadronic dynamical effects in the region of charmed hadron mass. 
This leads to overwhelming long-distance contributions and decreased reliability 
of theoretical predictions. Indeed, model-dependent evaluations of long-distance 
effects are possible (j6T |[62]) . which can be used to judge relative importance of 
long- and short-distance physics. 

These facts can constitute a problem for proper interpretation of new physics 
effects in FCNC processes. In addition, constraints on the strength of new interac- 
tions can be unambiguously placed only if the SM contributions are significantly 
smaller than the experimentally placed bound on a branching ratio. 

3.1.1 Inclusive and exclusive radiative decays c — >■ wy. Since rare 
radiative decays are two-body-decays, a branching ratio for exclusive or inclusive 
transitions is the primary observable. Thus, one has to evaluate relative NP/SM 
contribution for each model of NP. Only if the SM contribution, even dominated 
by the LD physics, is seen to be much smaller than current experimental bounds 
and possible NP contributions, such measurements can be useful in constraining 
NP models. Current theoretical estimates put decays rates of D° — > pj at the 
level of (0.1 -=- 0.5) x 10~ 5 and D° -> ^ at (0.1 -=- 3.4) x 10~ 5 . Currently, the 
decay D° — > (jyy has been measured to be (2.6to'g?io 17) x 10 5 (ISSJ) , and experi- 
mental constraints on other radiative decays are of the order of 10 -4 ()33|) . As the 
experimental bounds for radiative decays are pushed towards the SM theoretical 
estimates, these decays become less and less suitable to provide unambiguous 
constraints on New Physics models (|62p . 

In the SM, the radiative charm decays occur via the operators of the type 
O7 = (e/l6ir 2 )m c (ua^ u Pfic)F 11 " . In total, renormalization group running of per- 
turbative QCD requires a complete set of ten operators to describe this transi- 
tion (|6T]). Note that due to the chiral structure of the SM, the contribution of a 
similar operator O r 7 = (e/l6ir 2 )m u (ua^ u PLc)F llL ' is suppressed by a small factor 
m u /m c . Such suppression is not universal and is in fact absent in some models of 
NP, including SUSY. Thus, measurement of polarization of the final state photon 
can in principle be a nice probe of NP. 

3.1.2 Rare Decays D — > X u £ + £~. Decays of the type c —* u£ + £~ may 
allow a better separation of SM and NP effects. The simplest possible decay 
that is generated by this current is D° — > £ + £~ . Decays of this type are helicity- 
suppressed, with decay rates proportional to the masses squared of the final 
state leptons. This makes decays D° — ► e + e~ prohibitively small. Even the 
decay D° — > p, + p.~ is quite small. A calculation of short distance SM effects 
predicts a branching fraction of about 10 -18 (62). Long distance contributions 
bring the predicted branching fraction to an excess of 10 -13 , more precisely 2.6 x 
10~ 5 B(D° — > 77) (62). Thus, decays of this type provide almost background-free 
constraints on NP models. For example, R-parity-violating SUSY contributions 
are predicted at the level of 3 x 10~ 6 for some region of SUSY parameter space. 

In that sense, three body decays are more suitable for experimental studies, 
as they do not receive the above-mentioned helicity suppression. The two modes 
that have been studied most extensively are D — ► tt£ + £~ and D — ► p£ + £~ . These 
more complex final states provide additional tools to disentangle SM short and 
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long distance effects and NP phenomena. Figure 0] from Ref. (|62p illustrates 
this point with reference to the decay D + — * 7r + e + e~. It shows the predicted 
dilepton mass distribution normalized to T D +. The solid line represents the 
total SM prediction, while the curves represent predictions from a variety of 
minimal supersymmetric models. It is clear that for dilepton masses close to 
vector meson resonances such as p or there is no sensitivity to new physics 
contributions, however there are regions where NP effects are unambiguous. In 
particular, the region of low M e + e - is of great interest. Similar considerations 
apply to D — ► p£ + £~, where additional information is provided by the lepton 
forward-backward asymmetry 

„ f} dT / {dq 2 dx)dx — J , dT j [dq 2 dx)dx 

ApB{q } 55 dfJdJ 2 (13) 

where x = cos 6. Here 9 is the angle between the £ + and the D meson in the D 
rest frame. In the SM Afb(q 2 ) is negligibly small for all values of q 2 . 



3.2 Experimental Information 

A rare D radiative decay has recently been seen by Belle (65) with the reported 
branching ratio B(D° -> c/yy) = fifitolitofy x 1(T 5 . This branching fraction 
is measured by studying simultaneously the decays D° — ► (fry, D° — > <fiir°, and 
D° — > (jyq, as the last two modes, with higher branching fractions, induce some 
peaking background. D° — ► <jry is dominated by long distance effects, and the 
branching fraction is at the level expected from theoretical estimates. 

The most stringent limits on the di-lepton channels have been obtained recently 
by CDF fl66j), who study D° -> n+u~ and BaBar fl67|) who study both D° -» 
and D° — ► e + e _ . Both experiments use a D* tagged sample and normalize 
their results with respect to D° — > 7r + 7r _ . 

The channel c -> has been studied by CLEO-cflMl), BaBar (J67D, and DO 

(|68|) . The CLEO-c study focuses on L> + — > 7r + e + e~ because of their excellent 
sensitivity to di-electron final states, while BaBar studies both di-muon and di- 
electron final states of D + , Df and A+. DO has very recently reported results 
on the D + — > ir + fi + [i~ final state. All three experiments start by measuring 
D + — > 7r + — > -£ + £ _ 7T + to establish the order of magnitude of long distance 
effects. CLEO-c finds two events with an expected background of 0.02 events, 
BaBar finds 19 events over a background of 40 events, and DO finds 115 events 
over a background of 850 events. They find the branching fractions: 

B{D + -> e + e"7r + ) = (2.8 ± 1.9 ± 0.2) x 10" 6 [ (^4j) ] 

= (2.7±^)xl0- 6 [J5B], (14) 

and 

£(£>+ -> / u+^-7r+) = (1.8 ± 0.5 ± 0.6) x 10" 6 [dSHD] . (15) 

These experiments establish also 90 % confidence level upper limits on the short 
distance components of these branching fractions. Their results, compared with a 
representative sample of theoretical expectations based on NP scenarios are sum- 
marized in Table [6l These data show that experiments are reaching a sufficient 
sensitivity to limit the NP parameter space. The availability of higher statistics 
data samples from e + e~ experiments and collider experiments, as LHCb, when 
data taking starts, bear the promise of more stringent tests in the near future. 
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4 HADRONIC DECAYS 



Hadronic decays are interesting for several reasons. Absolute measurements of 
D meson branching fractions affect our knowledge of several D and B meson 
decays, from which CKM parameters are extracted. Multi-body final states 
provide information on light quark spectroscopy. Nonleptonic decays of charmed 
hadrons provide needed information (j72j) which helps in determinations of CKM 
angles (5 (173p and 7 (|74h in B decays and, can also help in determination of 
mixing parameters free from knowledge of hadronic strong phases (|75p . 



4.1 Theoretical Considerations 

Theoretical description of fully hadronic decays is significantly more complicated 
than leptonic or semileptonic ones, even though relevant effective Hamiltonians 
look quite similar to Eq. ([1]). Charmed nonleptonic decays are usually classified 
by the degree of CKM suppression. Least suppressed, where the quark level 
transitions are c — ► sud are labeled "Cabibbo favored" (CF) decays and governed 
by 

Hcf = ^V ud V c * s [Ci^)0 1 + C 2 ( f i)0 2 } + h.c, 

Ox = (siT^iukT^dk), 2 = {s i T^c k ){u k T^d l ) (16) 

where C n (/j,) are the Wilson coefficients obtained by perturbative QCD running 
from M\y scale to the scale \x relevant for hadronic decay, and the Latin indices 
denote quark color. 

The "Cabibbo suppressed" (CS) transitions are driven by c — > dud or c — > sus 
quark processes. Due to the presence of the quark-antiquark pair of the same 
flavor in the final state, the effective Hamiltonian takes much more elaborate 
form, 

H cs = ^^V uq v: q [C 1 ( f ,)Ol + C 2 ( f ,)Ol] 

q=s,d 

- ^K 6 y c * 6 ^<7 n ( M )0 + h.c, 

* n=3 

d = fer^Ci) (u k T>*q k ) , o 2 = (q i r IM c k )(u k r» qi ), 

where q = d, s, and O3-6 are the so-called "penguin" operators of the type 
{uc) V -A^2 q (qq)v±A (see, e.g. Ref. $76])). 

The "Doubly Cabibbo suppressed" (DCS) decay is the one in which c — ► dus 
quark transition drives the decay. The effective Hamiltonian for DCS decay can 
be obtained from Eq. (|16p by interchanging s «-> d. 

Calculations of hadronic decay rates governed by these transitions are quite 
complicated and model-dependent. Most often, simplified assumptions, such as 
factorization ([77][78|) are used to estimate the needed branching ratios. Some 
dynamical approaches, such as QCD sum rules, have been used to justify those as- 
sumptions (|79p . Charmed mesons populate the energy range where non-perturbative 
quark dynamics is active. This leads to resonance effects that affect the phases 
of nonleptonic decay amplitudes (|80p . Finally, standard methods of flavor SU(3) 
can be used in studies of non- leptonic D-meson decays (|8T|) . 
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4.2 Branching Fraction Measurements 

CLEO-c use tagged samples to obtain precise values for absolute hadronic branch- 
ing fractions for D° and D + (|82p , and for D s f)83|) . Their "double-tag technique" is 
similar to the one developed by Mark III ([Hi]) . From data at the -0(3770) they use 
three D° decay modes (D° -> R--K+, D° -» K-yr+vr , L>° -> 1^-^+7^+), and 6 
L>+ modes (L>+ -> i^-vr + vr + , L>+ -► i^-vr+vr+vr , L>+ -> K 5 7r+, D + -> ^vr+vr , 
_D + — > K57r + 7r _ 7r + , D + — > K + K~tt + . Single and double tag yields are used to 
extract the branching fractions and DD yields from a combined fit (85) to all 
the measured yields. This powerful technique, combined with careful efficiency 
studies based on data, dominates the present world averages (33). Corrections for 
final state radiation are included in these branching fractions. The large numbers 
of DD pairs, = (1-031 ± 0.017) x 10 6 and N D+D - = (0.819 ±0.012) x 10 6 

assure measurements at the 3% level, limited by systematic uncertainties. They 
apply a similar technique to derive absolute D s branching fractions from data at 
a center-of-mass energy of ~ 4170 MeV ()83|) . Here the dominant final state is 
D S D*, thus the analysis is more complex because of the 7 from the D*. 

B factories use D samples produced inclusively in B meson decays, or partially 
reconstructed samples of D or D s recoiling against a fully reconstructed charmed 
meson as normalization. Absolute branching fraction measurements for hadronic 
decays of D°, D + , and D s are summarized in Table which include new abso- 
lute BF's for D° -> K~tt + , reported by BABAR (86) with approximately a 2% 
uncertainty, and from Belle (I87p who measured the corresponding quantity for 
D s — > K + K~tt + with a precision of ~ 14%. 

In the BABAR measurement (j86j, B° -> D*+X£-i>i with D* + -» D°irf are 
identified by partial reconstruction. Events with a lepton i~ and a slow pion 
7T+ that could come from a D* + are selected by studying the reconstructed v 
invariant mass squared M^ meas inferred from conservation of energy and mo- 
mentum in the center-of-mass (CMS) system, using the D* 4-momentum inferred 
from tt s and the measured charged lepton's 4-momentum. The B normalization 
is determined from the peak in M^ meas (centered near for real D° events). 
Backgrounds under this peak, coming from D** , other BB combinations and 
continuum, are estimated from a wrong-sign lepton sample. Uncertainties in 
these backgrounds and in charged track reconstruction and particle identification 
efficiencies, dominate the systematic errors. 

The Belle result (|HT|) uses partially reconstructed e + e~ — > -D* + -D~ L (2536) events 
where D* + — > D^. They study two partially reconstructed samples: in the first, 
the D s is not reconstructed, but the 1)^(2536) is fully reconstructed in its decay 
to D*°K~ and D*~ K$- The soft 7 from the recoiling D* + is also required. In the 
other normalization sample, the D s is fully reconstructed in its K + K~ir + decay 
mode and is combined with the 7 to form a D* + . A recoil K from the D s \ decay 
was also required but the 1)^(2536) was not reconstructed. The result, obtained 
from the ratio of the K + K~n + signal in the first mode to the D s \ signal in the 
second, together with the (well-known) D* + — ► D°tt + and D* + — ► D s j BF's, 
provided the required BF for D s — > K + K~ir + . 

Of special interest are the decays to K . CLEO-c (88]) has recently studied 
both the K^tt and K$t^ final states. The Kl are identified as a peak in the 
missing mass. Effects of quantum correlations from the coherent D pairs from 
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"0(3770) decay are carefully taken into account. They measured the asymmetries 

(D) = B(D ^ K s n) - B(D ^ K L n) 
{ ' B(D -> K s tt) + B{D -> K l tt) 

for D = D° and D = D + . As pointed out in Ref. ([89]) D° -> if°7r° involves inter- 
ference between CF and DCS modes and, since we observe K° without knowing 
its strangeness, the K$ and Kl are related to give an asymmetry R(D) = 2 tan 9c 
where 9c is the Cabibbo angle. Rosner has observed (]90h that SU(3) flavor sym- 
metry, specifically [/-spin symmetry, predicts that the ratio of amplitudes for 
D° if °7T to D° -> K°tt° is tan 2 9 C ~ 0.054. This leads to the prediction 
that the value for D(D°) should be 2tan 2 f9 c ~ 0.109 ± 0.001. Ref. (JB8) finds 
D(L> ) = 0.108 ± 0.025 ± 0.024, significantly different from zero, and in good 
agreement with this prediction. There are no predictions for R(D + ), measured 
to be 0.022 ± 0.016 ± 0.018, compatible with zero. 

4.2.1 Cabibbo Suppressed Hadronic Decays. Due to CKM suppres- 
sion, these rates are expected to be lower by a factor r cs = | (V CS V US ) / (\V cs V u d) | 2 ~ 
0.05 relative to CF rates. 

Using their 281 pb _1 -0(3770) sample, CLEO-c (|9ip measured branching frac- 
tions for many multi-pion, r\ and to decay modes of D° and D + mesons. They use 
single tags and they extract absolute branching fractions using the correspond- 
ing well measured CF modes for normalization. These branching fractions range 
from (1 — 4) x 10 -3 and are measured with a precision of about (5 — 10)%. The 
largest rates are D° -> tt~tt + tt° (13.2 ± 0.6 x 10~ 3 ) for D° and D+ -> 7r"vr + vr + 
(3.35 ± 0.22 x 10" 3 ) and D + -> 7]ir+ (3.61 ± 0.36 x 10" 3 ) for D+. These mea- 
surements represent a significant improvement on previous knowledge, frequently 
being first observations. The results are generally consistent with simple CKM 
suppression. However, only an upper limit is extracted for D° — » 3-7T , in spite 
of the large B observed for D° — ► 7r~7r + 7r°. A possible explanation is that the 3 
pions are produced predominantly in an / = state, inaccessible to this mode, 
after a AI = 1/2 transition. 

The di-pion modes, D + — ► 7r + 7r , D° — > it + it~ and D° — > 7r°7r° are related 
by two amplitudes Aq and A 2 corresponding, respectively, to the 5-wave di-pion 
isospin 1 = and 1 = 2 states produced 

A+ ° = vl- 42 A+ ' = ^ + ^ 

A m = JlAt-^M (17) 

Following the procedure outlined in Ref. (|92p . CLEO obtains (|91|) from these new 
results \A 2 /Ao\ = 0.420±0.014±0.01 and avg(A 2 /A ) = (86.4±2.8±3.3)°, which 
is rather large. 

A long standing puzzle is found in the ratio 1Z(D° —* K + K~)/1Z(D° —* 
7r + 7r~) = 3.53 ± 0.12 where 7£ are BF's corrected for the 2-body phase space 
factor. 

BABAR (93]) has measured the branching fractions 

B{D° vr+Tr^vr ) = (1.493 ± 0.008 ± 0.055) x 10~ 2 
i3(L»° -> K+K--K ) = (0.334 ± 0.004 ± 0.015) x 10~ 2 . 
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The corresponding ratio for decays where an extra ir° is produced gives A4(D° — > 
K + K-ir°)/M(D° -> vr+vr-vr ) = 0.678 ± 0.027, in clear contrast to the 2-body 
ratio above. This result has recently been confirmed by Belle (|94p . 

CLEO-c also studied D s — > PP modes (f95|) . where "P" is any pseudo-scalar 
meson, using their = 4170 MeV sample. The modes (K + 7], K + n' , n + Ks, 
K + ir° are seen for the first time and are compared with their CF counterparts 
n + r], n + r]' and K + Kg. The ratios observed are reasonably consistent with the 
value of r cs « 5%. The decays £> s — ► 7r + 7r° have not yet been seen. The di-pions 
would be in an S-wave with 1 = 2, and would have to be reached through a 
AI = 2 transition, apparently much suppressed. 

4.2.2 Double Cabibbo Suppressed Hadronic Decays (c —* dus) of D 
Mesons. These decays are expected to be suppressed relative to CF modes 
by a factor r DCS = | (V c dV us ) / (\V cs V u d) | 2 ~ 3.1 x 10~ 3 . For Z?°'s, these rates 
are comparable to the mixing rate so that the two processes interfere; therefore 
disentangling the two effects needs some care. This is discussed in more detail in 
Section 

With such small branching ratios, one might wonder if DCS transitions could 
be affected by NP effects. However, since the final state is composed of quarks 
of different flavors, it is hard to find a well- motivated NP model that can affect 
DCS transition at an appreciable level (f96l) . 

The D° DCS branching fractions measured can be related, using world aver- 
ages (|33l) . to the expectations based on the value of r DCS 

(K+tt-) (1.45 ± 0.04) x 10~ 4 (1.18 ± 0.26) x r DCS 
(i^+vr-vr ) (2.96 ± 0.19) x 10" 4 (2.50 ± 0.57) x r DCS 
(K+tt-) (2.491°;^) x 10" 4 (2.10 ± 0.49) x r DCS . 

The only decay of D + , free from mixing effects, so far observed is D + — > K + tt°. 
BABAR dl7j) obtain B = (2.52±0.47±0.25±0.08) x 10" 4 and CLEO-c Qgfy obtain 
B = (2.28±0.36±0.15±0.08) x 10" 4 . Both experiments use the D+ -> K-tt + tt + 
as normalization, and the fourth uncertainty in each case is due to this. 

The decay of D + — » K + ir° has been observed by BABAR and confirmed with 
a more precise measurement by CLEO-c. Combining these measurements with 
the known lifetimes for D° and D + and the D° — ► K~ir + BF's from ref. ()33|) . a 
ratio can be formed, 

r (d+ -> k + tt°) 

-}— — — ( = (2.44 ± 0.33) x 10- 3 = (0.79 ± 0.20) x r DCS 

which is clearly compatible with the expected DCS rate. Similarly, the ratio of 
the two DCS rates is 

r (d+ -» k+tt°) 

— )—= f = 0.66 ± 0.09. 

r(p° A+vr~) 

A naive spectator diagram analysis would predict a ratio of 0.5. The difference is 
probably due to final state interaction effects, although annihilation or exchange 
diagrams could also contribute. 

4.3 Three Body Decays 

Multi-body D meson decays are a very rich source of information on long range 
strong interaction effects because of the complex interference patterns between 
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intermediate resonances formed between hadrons in the final states. In addition, 
some of them are relevant for CPV measurements in B decays (|74p or D° — D 
mixing parameters $73]). A vast body of experimental data has been studied, from 
fixed target experiments, CLEO (from data sets both near the T(45) and at 3770 
MeV center-of-mass energies), BABAR and Belle. Huge data samples - of order 
10 6 events in many channels - from the B factories have led to the need to review 
models used to fit the Dalitz plot distributions. Many results, no doubt, await 
such review and are yet to be published. Of particular concern is the extent to 
which phase information, essential to the determination of the parameters above, 
depends on the models assumed for these fits. Efforts are underway to attempt 
less model-dependent approaches (9!)ll 1001). 

4.3.1 Formalism for Three-Body D and D s Decays Decays of D or 
D s to three hadrons ABC often proceed through quasi 2-body modes D —* A + f 
followed by / —* B + C, where / is an intermediate "isobar" state, as outlined 
in Ref. (|77|) . When A, B and C are pseudo-scalar hadrons, the Dalitz plot, in 
which the squared invariant mass of one hadron pair is plotted against the squared 
invariant mass of one of the other pairs, contains all the dynamical information. 
These Dalitz plots often show intricate interference patterns between multiple 
resonances that may be produced in an intermediate state. 

The decay amplitude could be constructed from a partial wave expansion in 
any one of the three possible channels / defined by the particle pair. Here we 
choose / = BC. Each wave would then be characterized by the spin (J = L for 
pseudo-scalar hadrons) and isospin I of f 

oo 

A-M = ^2^2M L (p, q )F L j(s) (18) 

I L=0 

where s and s' are squared invariant masses for BC (i.e. /) and for the AC 
channel, respectively. F^j is the partial wave decay amplitude for the system 
/. M is a tensor function appropriate for the conservation of total spin in the 
decay and depends on L and the momenta p and q of B and A, respectively, that 
are defined in the / rest frame. The density of points on the Dalitz plot is then 
proportional to \A\ 2 . 

The complex function describes the production and final state scatter- 

ing resulting in the observed system /. Two distance scales may be distinguished. 
In the first, the parent D decays weakly, and hadronization occurs making an in- 
termediate hadron state k that may differ from /. At longer range, re-scattering 
(e.g. KK — ► 7T7r) occurs to make the observed system /. Thus, dropping label L 
and I, Fl j(s) can be written 

Ff(s) = T fk (s)Q k (s) (19) 

in which Tfk(s) is the matrix that describes hadron-hadron scattering. When s 
is small, T can only include elastic scattering. In this regime, therefore, in the 
absence of scattering between k and the recoil hadron A, the Watson theorem 
(jlOip requiring that the phase of Ff(s) should have the same s-dependence as 
elastic scattering, should hold. However, Qk(s) is an unknown function describing 
the short-range effects and could well have an s-dependent phase when strong k—f 
scattering takes place. In this case, the Watson theorem would not hold. A recent 
X-matrix fit to the 1 = 1/2 K~tt + S-wave amplitude from D + — > K~tt + tt + data 
by the FOCUS collaboration (102) indicate that this may be so. 



18 



Artuso, Meadows, and Petrov 



4.3.2 Analysis Methods of Three-Body D and D s Decays In fit- 
ting Dalitz plots, analysts have used several assumptions for the form of Ff{s). 
Most often used so far is the Breit-Wigner "isobar model" in which F(s) is ap- 
proximated as a linear combination of resonant terms for each wave. Eq. (I18p 
is approximated by a finite sum that includes only terms for the resonant states 
r observed in the data, no matter in which channel they occur. The amplitude 
(I18p is then a sum of Breit-Wigner propagators 



A(s,s') = NR + ^Ar 



A _ iS r M L (p,q)G L (q)G L ( P ) 
j±r — c r e 2 ■ : — — ■ {zu 
— s + im r L [s, L) 

with complex coefficients c r e lSr and a constant term NR, often introduced to 
describe direct, non-resonant decay to the three hadrons A, B and C. The par- 
ent D and the resonance have form- factors Gl that depend on L. Here, m r 
is the mass, and T(s,L) = T r (m r / \fs)(p/p r ) 2L+1 [GL{p) /Gl{Pt)] 2 is the mass- 
dependent width of the resonance r. The form-factors Gi{q) and Gi{p) for, 
respectively, the parent and the resonance r, are usually assumed to take the 
Blatt-Weisskopf form (|103p . Fractions are defined, for each resonance r, as 

Jf\A r \ 2 dsds' 

fr ~ fJ\NR + Z r A r \ 2 dsds'- (21) 

The sum of fractions, so defined, is not required to be unity, since the interference 
terms, included in the denominator, are missing from the numerator. Though 
Eq. (|2ip is in standard use, a better definition would be desirable. In particular, 
it works only to define resonant fractions from the isobar model. 

Other models attempt to address the problems associated with the description 
of the S- waves where the identity of resonances is less well understood. Un- 
derstandably, these methods do not have a convenient definition for resonance 
fractions. One model introduces hadron scattering through a ET-matrix i Uflll lOoj). 
related to the K and T-matrices and Q- and F-vectors discussed above by 

T kf (s) = (I-ipK(s))^K if (s) 
Q f (s) = Kji(s)P k (s) 

F f (s) = (I-ipK(s))- f l k P k (s) (22) 

where p is a matrix of phase space factors, purely imaginary below threshold, for 
any of the channels included. This is real, guaranteeing the unitarity of T, and 
it contains poles and non-resonant terms obtained from global fits to available 
scattering data. For example, Anisovich and Sarantsev (|106p derived a if-matrix 
representation of scalar tttt resonances through a global fit of all the available 
scattering data from threshold to 1900 MeV. The fit to the Dalitz plot then finds 
parameters for the production vector, Pk(s), which is, in essence, a function that 
describes the dependence of production of / on s. Pros and cons are known for 
this model. Unlike the isobar model, whose description of broad resonances by 
Breit-Wigner is flawed, this does preserve the unitarity of the T-matrix. It does 
not, however, do this for F, since the P-vector is arbitrary. Analyticity of T is 
not guaranteed either, so its ability to describe poles is limited. 
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Less model-dependent methods are also used. In restricted regions of the Dalitz 
plot, the angular "moments" of / can be used to measure the S, P, etc. ampli- 
tudes. A method that works over the whole Dalitz plot was introduced by the 
E791 collaboration (1107P in the analysis of D + — > K^tt + tt + decays. F(s) for the 
K~tt + S-wave was parameterized by a set of complex quantities at discrete s 
values. These quantities were treated as free parameters in their fit. The method 
requires a reference phase that, in their fit, was defined by the isobar model de- 
scription of the other waves. This introduced some model-dependence in their 
result. 

Assumptions about the models used to describe Dalitz plots affects results 
for CKM phases and D° - D mixing parameters. Extending such model- 
independence into fits becomes more important, therefore, as statistical precision 
of such measurements improves. 

4.3.3 Summary of Charm Dalitz plot analyses Several experiments 
have analyzed D° — > Kstt + tt~ (|33p . CLEO (|108p was first and included 10 reso- 
nances in their fit: K s p°, K s uj, K s / (980), ^/ 2 (1270), K s / 2 (1270), K s h (1370), 
K*(892)vr+, if*(1430)~7r-, K*(1430) vr+, K*(1680)tt + , and the doubly Cabibbo 
suppressed mode i^*(892) + 7r _ . CLEO found a much smaller non-resonant con- 
tribution than did the earlier experiments. The source of these contributions has 
been attributed to broad resonances such as the ifg(1430). The residual small 
non-resonant component may be a signature of broad scalar resonances such as 
the k and a. More recently, to extract information on the CKM angle 7(^3), 
analyses have also been made by both Belle ([10911110]) and BABAR( jllip using 
samples two orders of magnitude larger. At this level of statistical precision, 
the inadequacy of a simple isobar model to describe the data is revealed. After 
amplitudes for a clear signal for p — uj mixing and for radial excitations for the p 
are added, the fit is poor. Addition of a second, probably unphysical, a isobar, at 
~ 1000 MeV/c 2 improves the fit yet, even then, a reasonable fit quality is barely 
achieved. Alternate isobar models lead to uncertainties of about 10° in 7(^3), 
but these are smaller than the statistical and other systematic uncertainties of 
about 20° currently obtainable. 

Decays to D° — ► 7r + 7r~7r° are also used by BABAR to measure 7 (1112p . The 
Dalitz plot, also fitted earlier by CLEO (|113j) . is found to have a negligibly 
small NR component, and to be dominated by pn in all three charge modes. 
The structure of the plot shows strong, destructive interference between these 
modes in a six- fold symmetry suggestive of a dominant 1 = 0. This is consistent 
with the observation, noted in section 14.2. 1} that D° — > 3tt° decays are strongly 
suppressed. 

The decay D+ -» Tr+Tr+vr" was studied by E687, E791(P33I), FOCUS and, 
more recently, CLEO (]115p . This new CLEO analysis uses the largest sample 
(~ 4, 000 events) so far. An earlier isobar model analysis by E791 had reported 
the need to add a <r(500) Breit-Wigner to the -k + tt~ S-wave in order to get 
an acceptable fit. FOCUS re-examined this decay using a lf-matrix model for 
the 7r + 7r _ S'-wave, including hadron scattering data from a number of earlier 
experiments. The fit was acceptable, but the question of whether or not there was 
a er(500) was not clear. CLEO tried various other parameterizations for <r(500). 
Following a suggestion by J. Oiler (|116p a simple pole of the form l/(rn.Q — s) 
where tuq = (0.47 — 0.22i) was used rather than the Breit-Wigner used by E791. 
A scalar term based on the linear sigma model (j!17j) was also used. Both these 
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approaches produced a slightly improved fit. This all suggests that a low mass a 
resonance is likely to play a role in this mode, but even larger data samples are 
needed to find its pole parameters. 

E791 found evidence for a broad Kir scalar resonance in D + — > K + ir + ir + 
(I118p . Their original fit needed to include a non-resonant component with a fit 
fraction in excess of 90%. The inclusion of an S-wave Kir resonance with mass 
797 ± 19 ± 43 MeV and width 410 ± 19 ± 43 MeV improved the fit considerably 
and reduced the non-resonant fit fraction to 13 ± 5.8 ± 4.4 %. A new fit to this 
channel by Focus using three times the sample size and a i^T-matrix obtained 
from Kir scattering data indicates that an / = 3/2 Kir component is probably 
also present, but is unable to address the existence of a n. 

The charged .S-wave K^ir systems could provide new information on the 
k(800) seen, so far, only in the neutral K~ir + system. If it is an I = 1/2 
scalar resonance, it could also appear, with similar mass and width, in charged 
Kir systems like this. This was a motivation for BaBar to study the decays 
D° — > K + K~ir° (|119p . An isobar model fit is able to provide a satisfactory, 
though ambiguous, description of the data. In each model tried, the if* + (890) 
resonance is included among the components and is used as the reference for 
other phases. An S-wave K + K~ resonance is required to obtain a satisfactory 
fit, but there is no distinction between ao(980) and /o(980) at low masses. The 
higher mass K + K~ system also requires some contribution, and either /2(1525) 
or an fg with a similar mass works well. 

Three models are compared for the Kir S- waves. The first is a linear combina- 
tion of Breit-Wigner terms for k(800), K* (1430) and an NR term. The k,^ mass 
and width are allowed to vary. Second is the model that describes the data from 
the LASS Kir scattering experiment (|120p . The third uses the results for this wave 
from the E791 model-independent D + — * K~ir + ir + analysis (|107p . The second 
model gives the best fit, with the LASS phases shifted by ~ —90°. The fit with 
k is poor and requires a mass (870 ±30) MeV/c 2 and width (150 ±20) MeV/c 2 , 
the latter differing considerably from the width of ~ 400 MeV/c 2 reported for the 
neutral state. Establishing the existence of a k state with such a width requires 
more sophisticated analysis to find its pole in the T matrix so, though poor, this 
result does not rule it out as a genuine resonance. 

As the (j) band is relatively far from other interfering vector resonances, a 
study of the angular moments in the region around it is undertaken in order 
to attempt to learn more about the underlying 5-wave. The K + K~ invariant 
mass distributions are plotted, weighting each event by factors \J2l + ir/APi(8) 
(for £ =0, 1 and 2), where 6 is the angle between the K~ and the ir° in the 
K + K~ rest frame, and P% is the Legendre polynomial function of order I. These 
distributions, corrected for efficiency and with background subtracted, are defined 
as the "moments" in each mass bin. Assuming that only S- and P-waves 
(L = and 1) contribute to the K + K~ system, these moments are used to 
extract them: 

X = \^l±^l ■ X 1 = V2\S\\P\cose S p; A 2 = y||P| 2 (23) 

for each mass. In Fig. [5l the resulting magnitudes |5| and |P|, corrected for the 
phase space in the Dalitz plot (length of the K + K~ mass strips) are plotted. 
|P|, shown in (b), follows the <fi line shape well, with no asymmetry and little 
background, up to about 1040 MeV/c 2 . 
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In Fig. 0(a), \S\ values, determined from this channel, are compared with 
similar measurements of \S\ from an earlier analysis (1114p of D° — » K + K~ Kg 
decays. In that system, a^(980) in the K + K$ system was compared with the 
S'-wave in the K + K~ system and found to agree well, suggesting that ao(980), 
as opposed to /o(980) was present in both charge states. It is possible that the 
excellent agreement seen here is also evidence that ao(980) is the main contributor 
to the K~K + S-wave in this decay too. A more convincing test would, however, 
be to analyze the D° — » K~ir + r] and D° — ► r]ir°Ks systems since the oq would 
then be more obvious in its r/7r + decay modes. 

The last topic that we discuss is the Dalitz plot analysis of the decay D s —* 
K + K~ir + . E687 (|12ip reported the first Dalitz plot analysis of this channel 
with a sample of ~ 300 events. They found significant scalar contributions from 
/o(980) (or ao (980)). CLEO studied the K + K~ invariant mass spectrum in this 
decay (83) and found a clear peak at the <f> mass, but also a broad component with 
kinematic properties indicative of a scalar. This component is rather important 
because the mode D s — > (ftn has been commonly used for D s decay normalization, 
and the presence of this scalar component introduces an additional uncertainty 
in this branching fraction as the scalar channel that is absorbed in the (f)7r signal 
depends upon the experimental cuts used. For example, CLEO found that an 
uncertainty of the order of 5% was introduced depending upon the cut choices. 
BABAR (122) recently reported a preliminary fit of the Dalitz plot structure of 
this decay with over 100, 000 events. They studied moments in the low mass 
K + K~ and K~ir + systems which showed clear evidence for an S-wave contribu- 
tion in the former, but none in the latter. An isobar model fit shows, indeed, a 
strong /o(980) component, thus confirming that the Ds — > <j)ir branching fraction 
is not a very wise choice of normalization channel. 

A lot of work is still ongoing on the experimental side, with further exploration 
of modes relevant for CPV studies, both in D and B decays, and on the theory 
side to identify tools that reduce the model dependence in the Dalitz plot analyses, 
in particular when broad scalar resonances are involved. The rich structure of 
the Dalitz plots is not only a unique and powerful asset in understanding heavy 
flavor decay dynamics but also a tool for disentangling the intricacies of non- 
perturbative strong interaction effects. 

5 CHARM MIXING AND CP- VIOLATION 

The phenomena of mixing and CPV in the charm sector were first discussed three 
decades ago (|123p but the smallness of these effects are such that experimental 
evidence is scarce. CPV has not yet been observed, with upper limits currently at 
about the 1% level. On the other hand, after years of experimental investigation 
(|33p . evidence for mixing has finally been seen in two kinds of time-dependent 
measurements. The BABAR collaboration (124) has reported a 3.9 standard 
deviation effect in "wrong-sign" (WS) decays of D° — > K + tt~ 0. The CDF 
collaboration has reported a 3.8 standard deviation effect (|147p . Also, the Belle 
collaboration (145j) has reported a 3.2 standard deviation effect arising from the 
observed difference in lifetimes for decays to CP even final states D° -> K + K~ 
and 7r + 7r - compared to the mixed CP state K~ir + . This has also been confirmed 
by the BABAR experiment (j!48H . In addition, useful information on the strong 



1 Unless otherwise noted, charge conjugate states are assumed throughout this paper. 
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phase 5 affecting D° — » K + ir mixing resuits is coming from CLEO-c studies that 



exploit the quantum coherence of the D°TP pair produced near threshold ([33]) . 

5.1 Charm mixing predictions in the Standard Model 

The mixing arises from |AC| = 2 interactions that generate off-diagonal terms in 

the mass matrix for D° and D mesons. The expansion of the off-diagonal terms 
in the neutral D mass matrix to second order in the weak interaction is 

(m- 1 -t\ - 1 ID°\ H W=WU 1 ^ (D°\Hl ACl=1 \n)(n\Hl AC ^\D^) 
V 2 ) 2l 2Md 1 w 1 2Mb ~* M D -E n + ie 

(24) 

where H^= 2 and H^ =1 are the effective |AC| = 2 and |AC| = 1 Hamilto- 
nians. The off-diagonal mass-matrix terms induce mass eigenstates D\ and D 2 
that are superpositions of the flavor eigenstates D° and D , 



Di =pD°±qD° , (25) 

2 

where \p\ 2 + \q\ 2 = 1. The key quantities in D° mixing are the mass and width 
differences, 

AM D = Mi - M 2 and Ar D = T x - T 2 , (26) 
or equivalently their dimensionless equivalents, 

- AMd a - Ar ° 

xv = — , and y D = , (27) 

where Td is the average width of the two neutral D meson mass eigenstates. 
Two quantities, y§ p and y' D , which are actually measured in most experimental 
determinations of AFd, are defined as 

yo P = (r+ - r_)/(r+ +r_) = y D cos^ - x D sin <p - Ap ro ^j , 

x' = x D cos5 Kn + y D sin8 K7T , (28) 
y' = UD cos 5 K tt - xb sin 5 K tt , 

where the transition rates T± pertain to decay into final states of definite CP, 
Aprod = (j^D° ~ N-jf) I (^N D o + NjjOj is the so-called production asymmetry of 

D° and D° (giving the relative weight of D° and D° in the sample) and 5k-k 
is the strong phase difference between the Cabibbo favored and double Cabibbo 
suppressed amplitudes (f80|) . The quantities A m and cj> account for the presence of 
CPV in D°-D° mixing, with A m being related to the q,p parameters of Eq. (|25p 
as A m = \q/p\ 2 — 1 and <p a CP-violating phase of M 2 ± (if one neglects direct 
CPV) ([TO . 

The charm quark system is rather unique from the theoretical point of view, 
as its mass places it somewhere on the border of heavy and light quark systems. 
This makes prediction of D — D -mixing parameters a challenging task. As was 
shown in (|126|) . in the SM, xd and yo are generated only at second order in 
SU(3)i? breaking, 

x D ,VD~ sin 2 B c x [SU{3) breaking] 2 , (29) 
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where 6c is the Cabibbo angle. Therefore, predicting the SM values of xd and 
yD depends crucially on estimating the size of SU(3)jr breaking. 

Theoretical predictions of xd and yo within the SM span several orders of 
magnitude. Roughly, there are two approaches, neither of which give very reliable 
results because m c is in some sense intermediate between heavy and light. The 
"inclusive" approach is based on the operator product expansion (OPE). In the 
m c 3> A limit, where A is a scale characteristic of the strong interactions, AM and 
Ar can be expanded in terms of matrix elements of local operators (1127012811129} 
I130|,ll3ip . Such calculations typically yield xd, yD < 10~ 3 . The use of the OPE 
relies on local quark-hadron duality, and on A/m c being small enough to allow a 
truncation of the series after the first few terms. The charm mass may not be large 
enough for these to be good approximations, especially for nonleptonic D decays. 
An observation of yn of order 10 -2 could be ascribed to a breakdown of the OPE 
or of duality, but such a large value of yo is certainly not a generic prediction 
of OPE analyses. The "exclusive" approach sums over intermediate hadronic 
states, which may be modeled or fitted to experimental data (|132pi33jll34p . 
Since there are cancellations between states within a given SU(3) multiplet, one 
needs to know the contribution of each state with high precision. However, the D 
meson is not light enough that its decays are dominated by a few final states. In 
absence of sufficiently precise data on many decay rates and on strong phases, one 
is forced to use some assumptions. It was shown that phase space effects alone 
provide enough SU(3)i? violation to induce xn,yD ~ 10~ 2 (1261 1135]) . Large 
effects in yry appear for decays close to D threshold, where an analytic expansion 
in SU(3)i? violation is no longer possible; in addition, a dispersion relation can 
be used to show that in this case xd would receive contributions of similar order 
of magnitude. 

5.2 New Physics contribution to D° — D° mixing 

In order to see how NP might affect the mixing amplitude, it is instructive to 
consider off-diagonal terms in the neutral D mass matrix of Eq. (|24p . 

5.2.1 NP in |AC| = 2 interactions. Since all new physics particles are 
much heavier than the SM ones, the most natural place for NP to affect mixing 
amplitudes is in the |AC| = 2 piece, which corresponds to a local interaction 
at the charm quark mass scale. Integrating out NP degrees of freedom at some 
scale A, we are left with an effective Hamiltonian written in the form of series 
of operators of increasing dimension (ll36p . The complete basis of those effective 
operators, which most conveniently can be done in terms of left- and right-handed 
quark fields, is composed of eight operators, 



where Cj are the Wilson coefficients, and Qi are the effective operators, 




(30) 



Qi 

0,2 
0.3 



URJ^CRUL^CL 



Ob = URG^CLURa^CL 

Q& = ur-j^crur^cr, 
Q 7 = ulcrulcr, 



(31) 



ulcrurc l , 



URC L URC L , 



Q 8 = ULa^CRULCT^CR 
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Since these operators are generated at the scale (j, = A (at which NP is integrated 
out), a non-trivial operator mixing can occur if we take into account renormal- 
ization group running of these operators between fj, = A and [i ~ m c scales. This 
running can be accounted for by solving RG equations obeyed by the Wilson 
coefficient functions, 

-C( M )=7 T (aOC(m) 5 (32) 



dlog a. 

where 7 T (^) represents the matrix of anomalous dimensions of operators of 
Eq. (|3ip (136). A prediction for a mixing parameter x in a particular model 
of NP is then obtained by computing C(A) for a set of Qi(A) generated by a 
given model, running the RG equations of Eq.( [32]) and computing matrix ele- 
ments (D°\Qi(m c )\D ). 

Depending on the NP model, predictions for xn vary by orders of magnitude. 
It is interesting to note that some models require large signals in the charm system 
if mixing and FCNCs in the strange and beauty systems are to be small (e.g. the 
SUSY alignment model). A list of constraints on NP models is given in Table [8] 
(for more informative figures and methodology please see Ref. (|136p ). 
5.2.2 NP in |AC| = 1 interactions. The local |AC| = 2 interaction 
cannot, however, affect APd because it does not have an absorptive part. Thus, 
naively, NP cannot affect lifetime difference y. This is, however, not quite correct. 
Consider a D° decay amplitude which includes a small NP contribution, t4[D° — > 

n] = An M ^ + ^4n NP ^- Here, A„ P ^ is assumed to be smaller than the current 
experimental uncertainties on those decay rates. Then it is a good approximation 
to write y as 

VD - Ef'^^ + ^F 1 ^^ • (33) 

The SM contribution to y is known to vanish in the limit of exact flavor SU(3). 
Moreover, the first order correction is also absent, so the SM contribution arises 
only as a second order effect. Thus, those NP contributions which do not vanish 
in the flavor SU(3) limit must determine the lifetime difference there, even if 
their contributions are tiny in the individual decay amplitudes (|137p . A simple 
calculation reveals that NP contribution to y can be as large as several percent 
in R-parity-violating SUSY models or as small as 1CU 10 in the models with in- 
teractions mediated by charged Higgs particles (|137| PL38). 



5.3 CP-violation 

An observation of CPV in the current round of charm experiments is arguably 
one of the cleanest signals of physics beyond the SM (BSM). 

It can be easily seen why manifestation of NP interactions in the charm system 
is associated with the observation of (large) CPV. This is due to the fact that 
all quarks that build up the hadronic states in weak decays of charm mesons 
belong to the first two generations. Since 2x2 Cabibbo quark mixing matrix 
is real, no CPV is possible in the dominant tree-level diagrams which describe 
the decay amplitudes. CP-violating amplitudes can be introduced in the SM by 
including penguin or box operators induced by virtual 6-quarks. However, their 
contributions are strongly suppressed by the small combination of CKM matrix 
elements V c i,V* b . It is thus widely believed that the observation of (large) CPV in 
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charm decays or mixing would be an unambiguous sign for NP. This fact makes 
charm decays a valuable tool in searching for NP, since the statistics available in 
charm physics experiment is usually quite large. 

As with other flavor physics, CP-violating contributions in charm can be gen- 
erally classified by three different categories: 

(I) CP violation in the AC = 1 decay amplitudes. This type of CPV occurs 
when the absolute value of the decay amplitude for D to decay to a final 
state / (Af) is different from the one of corresponding CP-conjugated am- 
plitude ("direct CPV 11 ). This can happen if the decay amplitude can be 
broken into at least two parts associated with different weak and strong 
phases, 

A f = \Ax\ e iSl e** 1 + \A 2 \ e i&2 e^ 2 , (34) 

where fa represent weak phases (fa — > —<j)i under CP-transformation), and 
Si represents strong phases (Si — > Si under CP-transformation). This en- 
sures that CP-conjugated amplitude, Aj would differ from Af. 

(II) CPV in D° — D° mixing matrix. Introduction of AC = 2 transitions, either 
via SM or NP one-loop or tree-level NP amplitudes leads to non-diagonal 
entries in the D° — D mass matrix, 



r 

M-i- 



A p 



„ ^ *) (35> 

This type of CPV is manifest when R 2 m = \p/q\ 2 = (2M 12 - iT 12 )/(2M£ 2 - 

(III) CPV in the interference of decays with and without mixing. This type of 
CPV is possible for a subset of final states to which both D° and D° can 
decay. 



For a given final state /, CP-violating contributions can be summarized in the 
parameter 



X f = q -^ = R m e^ 
pA f 



A 



(36) 



where Af and Af are the amplitudes for D° — > / and D° — > / transitions 
respectively and S is the strong phase difference between Af and Af. Here <p 
represents the convention-independent weak phase difference between the ratio 
of decay amplitudes and the mixing matrix. 

Most of the experimental techniques that are sensitive to CPV make use of 
decay asymmetries, which are similar to the ones employed in B-physics f|139|, 

= T(D^f)-T(D^J) 

f r(D -> /) + t(d -> /)' 

One can also introduce a related asymmetry, 

r (D ^/)-r(^/) 
' r(D-/) + r(i>-/) 
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For charged .D-decays the only contribution to the asymmetry of Eq. (|37p comes 
from the multi-component structure of the AC = 1 decay amplitude of Eq. (|34"|) . 
In this case, 

2Im(A 1 A^)sm5 
Q/ | Ai| 2 + \A 2 \ 2 + 2ReA 1 A*cos5 

= 2rf sin 4>f sin 5, (39) 

where 5 = <5i — 62 is the CP-conserving phase difference and <fi is the CP- violating 
one. rt = \AzjAx\ is the ratio of amplitudes. Both ry and 5 are extremely difficult 
to compute reliably in D-decays. However, the task can be significantly simplified 
if one only concentrates on detection of NP in CP-violating asymmetries in the 
current round of experiments (|142p . i.e. at the 0(1%) level. This is the level at 
which (if is currently probed experimentally. As follows from Eq. (|39p . hi this 
case one should expect rj ~ 0.01. 

It is easy to see that the SM asymmetries are safely below this estimate. 
First, Cabibbo-favored (Af ~ A ) and doubly Cabibbo-suppressed (Af ~ A 2 ) 
decay modes proceed via amplitudes that share the same weak phase, so no 
CP-asymmetry is generated!]. Moreover, presence of NP amplitudes does not 
significantly change this conclusion (|96H . On the other hand, singly-Cabibbo- 
suppressed decays (Af ~ A 1 ) readily have two-component structure, receiving 
contributions from both tree and penguin amplitudes. In this case the same 
conclusion follows from the consideration of the charm CKM unitarity, 

v ud v; d + v us v* + v ub v; b = o. (40) 

In the Wolfenstein parametrization of CKM, the first two terms in this equa- 
tion are of the order 0(A) (where A ~ 0.22), while the last one is 0(A 5 ). 
Thus, CP-violating asymmetry is expected to be at most af ~ 10 -3 in the 
SM. Model-dependent estimates of this asymmetry exist and are consistent with 
this estimate (|143p . Other observables are also possible, such as untagged CP- 
asymmetry (|144p . which for some final states can be written in terms of experimentally- 
measured quantities only, i.e. which suffers no theoretical uncertainties. 

5.4 Experimental Considerations 

A number of experimental realities prevented positive evidence for D°D° mixing 
until last year when sample sizes, with sufficiently small backgrounds finally be- 
came available (|1241ll45j) . Among the most important was that the mixing rate 
was very small (R M = (xl + y 2 )/2 ~ lO" 3 ) and mixing effects were greatest at 
about two lifetimes where sample sizes were quite depleted. Backgrounds were 
high and use of the two way constraint on both invariant mass M of the decay 
products of the D° and the difference, AM between M and the D* invariant 
mass used to tag the D° flavor was necessary to suppress these. Cuts in data 
samples aimed at removing D°'s from B decays that adversely affect decay time 
distributions also ate into sample sizes. Also, the D° lifetime (~ 400 fs) was small 
on the scale of measured decay times t and their uncertainties 5t and large sam- 
ples were required to properly understand their resolutions and the time offsets 
arising from detector alignment issues. 

2 Technically, there is a small, C(A 4 ) phase difference between the dominant tree T amplitude 
and exchange E amplitudes. 
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5.5 Experimental Studies of D° — D° Mixing 

Studies of D° - W mixing mostly use samples of D° mesons whose flavor is 
identified by the sign of the "slow pion" tt s in the decay of a D* + — > D°Trf. 
Effects of mixing in the subsequent decays are then examined in one of two 
alternative ways. Studies of time-dependent decay rates look either for structure 
in the "wrong sign" (WS) hadronic final states or for differences between the 
lifetimes of decays to CP eigenstates and to states of mixed CP. Alternatively, 
time-integrated rates are studied looking for decays either to semi-leptonic final 
states with the WS lepton, H or for the effects upon rates resulting from quantum 
correlations between D° — D systems produced in ifi(2S) decays at CLEOc. 

5.6 Time-Dependent Studies 

Mesons produced at time t = as D°(D°) have amplitudes for decay to final 
state /(/) at time t given by: 



(f\H\D°(t)) = e - {TD+iMD)t A f 
(f\H\lf{t)} = e - {rD+iMo)t Af 



cosh((y D + ix D )—) + X f sinh((y D + ix D )-^-) 

r t r t 

XJ 1 sinh((y D + ix D )-^-) + cosh((y D + ix D )-^-) 



(41) 



where Xf = is defined in Eq. (|36p . M D and T D are the mean mass and width 

of the D\ and D2 and Af and Af are amplitudes describing, respectively, the 

direct decays of D° and D to final state /. 

In the absence of CPV in either mixing or direct decay to /, then Xf = e lS 

where 5 is the strong phase difference between D° — > / and D — > / amplitudes. 
With CPV in mixing (p ^ q) or in direct decay (\Af/Af \ ^ 1), Xf will be complex 
with weak phase cpf and a strong phase 5. 

Decay rates, proportional to the square modulus of these amplitudes provide 
information on x D and y u , only if the strong phase 5 is known. This is, however, 
zero in the important special case when / is a CP eigenstate. In that case, if / 
has CP = rj the direct decay amplitudes are related by A v = rjA v , so that 5 = 0. 

5.6.1 Results from D° — > K + tt~ Decays. These WS decays can take 
place in two ways - either directly, by a doubly Cabibbo suppressed (DCS) mecha- 
nism or by mixing (D° -> D ) followed by "right sign" (RS) Cabibbo favored (CF) 

decay D — > K + tt~ . These two processes interfere giving the time-dependent de- 
cay amplitudes given in Eqs. (|4ip . Neglecting CPV and assuming that \x D \ and 
\y D \ <C 1), this leads to the decay rate Rws for these decays 

oc R D + ^R~y't + R M {T u tf/2 (42) 

R D is the DCS decay rate alone, in the absence of mixing, R M is the mixing rate 
and the middle term, linear in y', results from interference between mixing and 
DCS amplitudes. Deviation from a purely exponential decay, expressed in the 



3 Actually, the signal for WS decays also do use information from the time-dependence from 
mixed decays to distinguish them the backgrounds. 
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right side of Eq. (|42p . can provide information on x' 2 and y', but not the values 
of x D and y D or their relative signs. 

This method has been used a number of times, in E691 and E791 (fixed target 
experiments at Fermilab), by BABAR and by CLEO using 9 fb _1 of e + e~ col- 
lisions at yfs ~ 10 GeV (|146p . The most stringent limit on D°D mixing using 
this method is reported by Belle (J149J). Recently, however, using a 384 ftr 1 data 
sample, the BABAR collaboration finally reported evidence for mixing from their 
sample of 4, 030 ± 90 decays of D° mesons to the WS final state. 

The distribution in decay times t (uncertainty 5t) for these D° and D events 
is shown in Fig. EJa). Also shown are the projections onto the time axis of fits 
made in the (M, AM, t, 5t) distribution for events near the D° signal region. The 
dashed curve shows a fit made on the assumption that x' 2 = y' = and the 
solid curve represents a fit where these parameters are allowed to float. Precise 
knowledge of parameters describing the distribution of 5t, clearly an important 
component in these fits, is obtained from a sample of ~ 1.3 x 10 6 "right sign" (RS) 
candidates for the Cabibbo favored decays D° — > K~ir + . The difference between 
the fit with mixing and that without is significant as seen from the residuals in 
Fig. E](b) and in the ratio of WS to RS decays in Fig. [6^c). The central values 
from the mixing fit occur at a negative (unphysical) value for x' 2 . The likelihood 
obtained in the fit at the point (x' 2 = 0, y' = 6.4 x 10 -3 ) just inside the physical 
region differs from that at the central point by only 0.7 units, while the likelihood 
at the no mixing point {x 12 = y' = 0) differs by almost 24 units. This difference is 
taken to indicate evidence for mixing at 3.9'cr' significance. The CDF experiment, 
using data from a 1.5 fb" 1 pp exposure at yfs = 1.96 TeV, has recently confirmed 
the BABAR result at a significance of 3.8 standard deviations (|147p . These results 
are also summarized in Table [9l where it is seen that agreement with BABAR 
is excellent. In this experiment, it is not possible to remove background from 
D mesons from B decays simply with a momentum cut, as in the BABAR and 
Belle experiments. This is, however, identified from its vertex distribution and 
properly taken into account in their analysis. In an earlier analysis, the Belle 
collaboration (149) also obtained results with greater precision than either the 
BABAR or CDF measurements. Their central values for x' 2 and y', however, 
were closer to zero and so did not provide convincing evidence for mixing and 
only 95% confidence limits could be reported. 

Table [9] summarizes the data on mixing derived from this hadronic channel. 
Generally, the highest sensitivity is achieved with fits assuming CP violation. 
Note that the parameters x' and y' include a strong phase 5, the phase of the 
amplitude ratio (K + ir-\D°) / (K+ir-\D°). 

5.6.2 Lifetime Difference Measurements. The Belle collaboration also 
presented evidence for mixing using a 540 fb" 1 e+e" data at the T(45). They 
measured lifetimes for singly Cabibbo-suppressed decays to the CP = +1 final 
states K + K~ and 7r + 7r~ and for Cabibbo-favored decays to the final state K~ir + 
with mixed CP flUS]). 

This method was first used by the E791 collaboration (|150|) who were unable 
to detect mixing due to limited statistical precision. Subsequently, upper limits 
were reported by FOCUS, and CLEO. In the approximation that x D and y D are 
small, decays to CP = +1 eigenstates follow approximately exponential forms 
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with lifetimes (j!25j) . respectively, for decays of D° and D of 

t + = r° [1 + \q/p\(y cos (f)f - xsin^)] -1 

r~ = r° [1 + \p/q\(y cos <fif + x sin (43) 

where T° is the lifetime for decays to non-CP eigenstates such as K 7r + . For 
such measurements, it is convenient to define 

AY = (t°/(t>K (44) 

where (r) is the average of t + and r _ and a r = (r~ — t + )/(t~ + r + ) is their 
asymmetry. In the absence of mixing {x = y = 0) both are zero. In the absence 
of CPV in mixing or in decay (i.e.(f>f = 0) then Ay = and y£ p = y. 

These measurements require that backgrounds are small and have a well un- 
derstood time-dependence. The Belle samples, consisting of 111K K + K~, 1.22 x 
10 6 K~tt + and 4QK vr+vr- with purities 98% and 99% and 92% purity, respec- 
tively. The decay times for these were fit simultaneously to distributions with an 
exponential for each signal convolved with the time resolution function over the 
expected background distributions. 

Estimates for y^ p were made using both D° and D samples together and 
for t + and r~ from separate fits to each. The major systematic uncertainties 
were from an understanding of time offsets and from possible i-dependence of the 
efficiency for reconstructing events. The result obtained 

yZ P = (1.31 ±0.32 ±0.25)% 
a T = (0.01 ±0.30 ±0.15)% 

shows that y° p is not zero, evidence for mixing at the 3.2cr level. However, a T is 
consistent with zero, so there is no evidence for CPV. 

A similar analysis by BABAR, using a 384 fb _1 data set, confirms these results 
(|145p . Though the BABAR sample is smaller, a higher purity is achieved. Their 
results 

y° p = (1.24 ±0.39 ±0.13)% 
Ay = (-0.26 ±0.36 ±0.08)% 

agree well with Belle's, and show evidence for mixing at the 3. Oct level, but with 
no evidence for CPV. 

5.6.3 Mixing in the Decays D° — * K s tt + tt~. For D° —* K S TT + it~ , the 
time dependence of the Dalitz plot distribution allows one to measure x D and y D 
directly. This technique was first developed by the CLEO collaboration (151), 
who used a 9 fb _1 data sample. Assuming no CPV , they obtained the limits 
(-4.7 < x < 8.6)% and (-6.1 < y < 3.5)% at 95 % confidence level. 

A special feature of this decay mode is that, treating K$ as a CP = ±1 
eigenstate, final states / reached via the tttt channel are also CP eigenstates 
for which the strong phase difference 5 = 0. This provides a reference phase 
for a time-dependent Dalitz plot analysis that will determine the At, At for all 
channels everywhere in the plot. Thus, these decays allow measurement of x D , 
y D and their relative sign, free from any unknown strong phase. This fit can also 
provide magnitude and phase of the Ay's, providing a test for CPV. 

Belle applied this technique to a data sample 60 times larger. The Kstt + tt~ 
Dalitz plot is shown in Fig. [71 where the RS K*~ form a vertical band and 
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p° is the diagonal band. The WS K* + appears as a weak, horizontal band 
destructively interfering with other structures in the plot. It is expected that 
RS and WS K* amplitudes should have opposite signs since their weak CF or 
DCS phases differ in sign. An isobar model description of the amplitudes At 
is used, writing them as sums of Breit-Wigner functions and spin factors, with 
complex coefficients determined in the fit as in Eq. (I20p . For the Belle data, 18 
isobars, including p — to mixing, two ir + 7r~ 5- wave states and a non-resonant term 
NR were required to provide a reasonable match to all the features seen. The 
appropriate time-dependencies were included as in Eq. (|4ip . 

Three fits were made. In one, no CPV is included (all D° and D events 
were combined and isobar coefficients for D° set to equal those for D . The 
condition p = q is also imposed. In the other fit these conditions were relaxed, 
introducing a set of complex isobar coefficients for D differing from those for D°. 
The modulus and phase of the ratio p/q were also allowed to float. In the last 
fit, CPV in mixing (p ^ q) is allowed but direct CPV is not (isobar coefficients 

constrained to be the same for D° as for D ). The results of these three fits are 
summarized in Table fTUl 

The authors estimate that the best solution differs from the mixing point (x D = 
y D = 0) by 2.2 standard deviations. Allowing for CP violation, they obtain the 
CPV parameters \q/p\ = 0.86±g;|g±g;g| ± 0.08 and avg(q/p = (— 14±J|±|±t) - 
This result does not establish evidence for mixing nor for CPV. It does, however, 
illustrate a powerful way to study mixing, able to determine x D , y D and their 
relative signs, and the CPV parameters (f) and \p/q\- 



5.7 Time Independent Studies. 

Rates for WS leptons in semi-leptonic decays, a clear signal for mixing, could 
directly determine the mixing rate R m . So far, with samples available, these 
have only been able to produce upper limits on this quantity. The most precise 
limits are from Belle (152[), from a 253 fb _1 sample, and from BABAR (153) with 
a 344 fb^ 1 sample. The BABAR analysis differed considerably from Belle's in 
that they adopted a double tagging approach, requiring a fully reconstructed D 
on the side opposite the semi-leptonic decay, used to reduce WS backgrounds. 
Event yields differed by orders of magnitude between the two experiments, yet 
the limits obtained were very similar: 

Belle R M <lx lO" 3 

BABAR (-1.3 < R M < +1.2) x 10~ 3 

CLEO-c can exploit the quantum coherence of the D°D° produced via the 
V>(3770) resonance to extract several important variables affecting D°D mixing. 
The pair is produced in a C = — 1 state, while C = +1 states are accessible 
at a higher energy, if ^yD°D final states can be tagged. By reconstructing one 
neutral D meson into a CP eigenstate decay mode, the CP eigenvalue must be 
the opposite for C odd wave functions, and the same for C even wave functions 
(|154pi55|) . assuming no CPV. All the mixing observables can, in principle, be 
measured. A sophisticated fitting technique has been developed (|156p to reach 
maximum sensitivity. For decays of ^(3770) to D° — D° pairs (|157|J158|) . the D 0, s, 
in a coherent P-wave state, and with opposite flavor and CP, decay in a correlated 
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way so that K ir + rates depend upon R M , R n y' and 6, the strong phase difference 
between amplitudes for CF and DCS decays. Using their 281 pb _1 data sample 
collected at the ^(3770), CLEO determines cos 5 = 1.03±°;^ ± 0.06 (fl59l) . 

5.7.1 Averaging the Results. The decays discussed above provide infor- 
mation in different forms, depending upon the final state. Semi-leptonic modes 
determine R M , WS hadronic systems measure x' 2 and y' separately for D° and 
D . Decays to CP eigenstates measure y^ p and Ay and quantum correlated 
states from -0(3770) — ► D° decays can measure R M , R D , y D and cosc/> for various 
hadron systems. Parameters x D , y D , \p/q\ and avgp/q are obtained from time 
dependence amplitude analyses of decays of D° to final states with more than 
two hadrons, as long as those amplitudes include at least one CP eigenstate. 

The parameters of physics interest that define values for all these quantities 
include x D , y_ D , \p/q\, aigp/q, , <f%S , R D and its asymmetry a D = (R D - 
R D )/(R D + R D ). The Heavy Flavor Averaging Group (HFAG) made a y 2 fit to 
obtain values for these parameters that best describe all 26 available observations. 
(|160p Results, projected onto the (x D ,y D ) and (|p/<z|, argp/g) planes in Fig. [HJ 

The point at which there is no mixing on the (x D ,y D ) plane lies at the origin, 
outside the 5a contour. This indicates very strong evidence for D° — D mixing. 
There is no evidence, however, for CPV in the mixing. The point where p = q in 
the (\p/q\, aigp/q) plane lies at (1,0), right on the la contour. 

Evidence for mixing is convincing and different experiments and methods agree 
well. Since y£ p > 0, then the CP = — 1 state has longer lifetime, as in the 
other neutral mesons that mix. If the sign of x D /y D remains positive as more 
measurements are made, then the CP = —1 state is lighter, unlike the K° system. 
Finally, there is no evidence for CPV so far. 

5.7.2 Results of Experimental Searches for CPV. At the present 
time, there is no experimental evidence for CPV in weak decays in the charm 
sector. Finding it in CF or DCS decay modes would signal possible NP. For decays 
to multi-body systems, the decay asymmetry could appear only in certain sub 
channels and result in particle-antiparticle differences in phase space distribution. 
Searches may also be made for T-violation. 

5.7.3 Asymmetries in Time-integrated Partial Widths. Until this 
year, measurements of the asymmetries in partial widths for time- integrated D 
decay rates were known with precisions of a few percent. A recent compilation of 
average values f|161 1) for D decay asymmetries for a number of modes measured 
by CDF, FOCUS, E791, CLEO, BABAR and Belle are summarized in Table EU 
They are all consistent with zero, with a precision of a few % in most cases. 
Included in the most recent results are measurements from CLEOc ([6]) of CF 
modes using their -0(3770) data. Precisions were excellent, of order 1%. 

The most precise measurement comes from BABAR f!162j) with a precision of 
0.34%,for the decays D° -> K+K~ and of 0.5% for D° tt+tt". This pre- 
cision required careful consideration of systematic effects, notably charge and 
tagging asymmetries calibrated using data from the K~ir + CF mode. Effects 
from forward-backward production asymmetry arising from higher order QED 
were also taken into account. Systematic uncertainties were 0.13% for K + K~ 
and 0.22% for 7r + 7r~, so the results were limited by statistics. Nevertheless, find- 
ing CPV closer to 0.1% in CS modes, and less in CF or DCS modes (the goals 
to observe NP) will require even better precision. 

CPV effects would influence the Dalitz plot distributions in three-body decays 
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such as D° — > 7r + 7r 7r since asymmetries would affect different partial waves 
in each of the three channels in differing degrees. Such effects would also be 
expected to introduce phase differences between D° and D , and these could be 
observed in analysis of the Dalitz plots. ()113[) found no such effects at the few % 
level. A model-independent approach might be interesting to pursue - possibly 
by the B factories where an order of magnitude more data is available. An isobar 
model analysis of this mode by CLEOc sees no effect at the few percent level. 

5.7.4 T- violation Studies. It has been pointed out (I139P that T-violation 
in charmed meson decays may also be observed in the asymmetry of triple scalar 
products (T-odd) of the momenta of the particles emerging from 4-body de- 
cays (139). The only such measurements to date are from the FOCUS collabo- 
ration (I163h . who measured the quantities (defined in terms of momenta p with 
suffices indicating each product from 4-body D°, D and Df decays) 

Ct = PK+ ■ (Ptt+ x Ptt-) 

T(C T > 0) - T(C T < 0) 
T T(C T > + T(C T < 0) 

The conjugate quantities (Ct and At) for D° , D~ and Dj decays were also 
measured. The modes studied were D — ► K + K~tt + it~ , D + — ► K + Ksir + ir~ 
and — » K + Kstt + tt~ . Any asymmetry not consistent with zero would indi- 
cate T-violation in the absence of strong interactions. The latter introduce the 
same asymmetries in particle and anti-particle decay, so the quantity ^T-viol = 
(1/2) (At — At), which would be zero if T is conserved, is evaluated 

D° K+K-tt+tt- 0.010 ±0.057 ±0.037 
D+ -> K + K s ir + ir- 0.023 ± 0.062 ± 0.022 
D° K+K-tt+tt- -0.036 ±0.067 ±0.023 

all consistent with zero. 

5.7.5 Summary. Clearly, experimental precision is not yet sufficient to chal- 
lenge the SM with respect to its predictions of CPV in the charm sector. The 
outlook is, however, good since the precision of asymmetry measurements is not 
yet limited by systematics. Also, model-independent studies of the multi-body 
channels 7r + 7r~7r and K + K~~tt°, less prone to systematic uncertainties, show 
promise as a tool for observing effects of CPV. 

Both BABAR and Belle still have large samples of 4-body decays where T- 
violation tests similar to that made by FOCUS can be repeated. Also, more data 
is to come from both B factories on the K~K + and 7r~7r + channels where the 
precisions are beginning to become interesting. Even larger samples from charm 
factories, LHCb, or, possibly, Super B factories in Italy or at KEK might produce 
a definitive answer on NP manifest in these studies. 

6 CONCLUSIONS AND OUTLOOK 

Charm decays remain an exciting field for both theoretical and experimental 
investigations. Charm quark transition amplitudes, described in this review, 
represent a crucial tool to understand strong interaction dynamics in the non- 
perturbative regime. Complementary information that constrains model building 
and lattice gauge calculations is coming from the rich spectroscopy of charmed 
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mesons and baryons, which is beyond the scope of this review. The validation 
of theoretical tools that tackle non-perturbative processes is critical to precision 
tests of the Yukawa sector of the SM, in particular to unitarity checks of the 
Cabibbo-Kobayashi-Maskawa matrix. 

Finally, charm decays provide a unique window on NP, provided it affects 
u-type quark dynamics. This way charm phenomenology can have an impact 
on the interpretation of results from the direct searches for new physics to be 
performed at the LHC. Charm quark is the only -u-type quark that can have 
flavor oscillations. Thus the observation of D°D mixing is already constraining 
many scenarios of physics beyond the SM. In addition, a multitude of new physics 
models predict enhancements on CP violating phases in D decays, beyond the 
10~ 3 level generally predicted within the SM. A full exploration of CPV in the 
charm sector, that hopefully will be achieved in the next decade, is a critical 
ingredient to further narrow the vast parameter space presently characterizing 
all the new physics models. 
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Table 1: Theoretical predictions for f D +J D+ J D +/f D+ . 

Authors f D+ (MeV) / fl +(MeV)' f D +/fp+ 

Unquenched lattice calculations 

HPQCD+UKQCD JHJ 208 ±4 241 ± 3 1.162 ±0.009 

FNAL+MILC+HPQCD ([TOD 201 ± 3 ± 17 249 ± 3 ± 16 1.24 ± 0.01 ± 0.07 

Quenched Lattice QCD Calculations 

Taiwan (20) 235 ± 8 ± 14 266 ± 10 ± 18 1.13 ± 0.03 ± 0.05 



UKQCD (j2U) 210±10l^ 236 ±St\l 1.13 ± 0.02+°;°^ 

Becirevic et al.fl22]) 211 ± 14±^ 2 231 ± 12+^ 1.10 ±0.02 

QCD sum rules and other approximations 



J. Bordes et al.mh 177 ± 21 205 ± 22 1.16 ± 0.02 ± 0.03 

S. Narison dMI 203 ± 10 235 ± 24 1.15 ±0.04 

Field Correlators (25]) 210 ± 10 260 ± 10 1.24 ± 0.03 

Isospin Splitting (j26} 262 ± 29 



120 




0.50 



MM 2 (GeV 2 ) 



Figure 1: CLEO-c MM 2 using D~ tags and one opposite charged track with 
no extra energetic clusters. 10 The insert shows the signal region for D + — > av^ 
enlarged; the defined signal region is shown between the two arrows. 
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Table 2: Results for B^ = B(D S — > £>(-D s — > T + z/ r ); and / D +. (Numbers 

have been updated using L> s lifetime of 0.50 ps.) Results below the line have not 
been used in this average from Ref. (115p . The assumed value of B(Df — > (j)ir + ) 
is listed whenever available. ALEPH average their two results to obtain a value 

for fp+- 



Exp. Mode B(xW*) f D + (MeV) 



CLEO-c 


p+v^ (34) 


5.94 ±0.66 ±0.31 




264 ± 15 ± 7 


CLEO-c 


T+U T ([MI) 


80.0 ± 13.0 ±4.0) 




310 ± 25 ± 8 


CLEO-c 


r+z/ r (J35D 


61.7 ±7.1 ±3.6 




275 ± 10 ± 5 


CLEO-c 


combined 






274 ± 10 ± 5 


Belle 


P + v^ 


6.44 ±0.76 ±0.52 




279 ± 16 ± 12 


Average 








275 ± 10 


CLEO dHTJ) 




6.2 ±0.8 ±1.3 ±1.6 


3.6 ±0.9 


273 ± 19 ± 27 ± 33 


BEATRICE ([35 




8. 3 ±2. 3 ±0.6 ±2.1 


3.6 ±0.9 


312 ± 43 ± 12 ±39 


ALEPH ([39D 




6.8 ±1.1 ±1.8 


3.6 ±0.9 


282 ± 19 ± 40 


ALEPH dM]) 


T + V T 


58 ± 8 ± 18 






L3 (@0D 


T + U T 


74 ± 28 ± 16 ± 18 




299 ± 57 ± 32 ± 37 


OPAL (Hp) 


T + V T 


70 ± 21 ± 20 




283 ± 44 ± 41 


BaBar (28) 




6.74 ±0.83 ±0.26 ±0.66 


4.71 ± 0.46 


283 ± 17±7± 14 



Table 3: Summary of recent absolute branching fraction measurements of exclu- 
sive D + and D° semileptonic decays. When only the CLEO-c absolute number 
is available, no average number is provided. 



Decay 


mode 


B(%) [CLEO-c] (31) 


B(%) [BES](32) 


B(%) [PDG06 average] 




K~e + u e 


3.44 ±0.10 ±0.10 


3.82 ±0.40 ±0.27 


3.47 ±0.13 


D ^ 


ir~e + v e 


0.262 ± 0.025 ± 0.008 


0.33 ±0.13 ±0.03 






K*~e + i/ e 


2.16 ±0.15 ±0.08 






D° — > 


p~e + v e 


0.194 ±0.039 ±0.013 








K"e + u e 


8.71 ±0.38 ±0.37 


8.95 ±1.59 ±0.67 


8.7 ±0.5 


D+ -» 


TT°e + U e 


0.44 ±0.06 ±0.0.03 






D+ -» 


K*°e+v e 


5.56 ±0.27 ±0.23 






D+ -. 


p°e + u e 


0.21 ±0.04 ±0.01 








■ we + v e 


0.16i°;$:±0.01 







Table 4: Fit parameters in Eq. (8), decay rates and CKM matrix elements. The 
first errors are statistical; the second systematic. 



p 


F 


a 


P 


7T 


0.64(3)(6) 


0.44(4) (7) 


1.41(6) (13) 


K 


0.73(3)(7) 


0.50(4)(7) 


1.31(7)(13) 
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Figure 2: Fits (solid lines) to the U distributions in CLEO-c 28 data (dots with 
error bars) for the five D + semileptonic modes: (a) D + — > K°e + i , e , (b)D + — > 
K*°e + f e; (c) L>+ -> vr°e + ^ e , (d)D+ -> p e + i/ e , (e)D+ we + i/ e . The arrows in 
(e) show the signal region. The background (in dashed lines) is visible only in 
(c) and (d). 



Table 5: Measured shape parameter a compared to lattice QCD predictions. 



a 


(D u - 






Lattice QCD 




0.5 ±0.04 ±0.07 


FOCUS (jSU) 




0.28 ±0.08 


±0.07 


CLEOIII ([Ml 




0.36 ±0.10 


-1-0.03 
-0.07 


Belle © 




0.52 ±0.08 


±0.06 


BaBar J53J) 




0.38 ±0.02 


±0.03 


a(D u - 






Lattice QCD 


(H2D 


0.44 ± 0.04 


±0.07 


CLEO III (|53 


J 


o 7 +0.20 ■ 
u -°'-0.31 ^ 


0.15 


Belle © 




0.10 ±0.21 


±0.10 
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Figure 3: Belle from factors for (a) D° — > K 



q 2 (GeV 2 /c 2 ) 



in q 2 bins of 0.067 GeV 2 



and (b) D —* ir l + u in q 2 bins of 0.3 GeV 2 . Overlaid are the predictions of 
the simple pole model using the physical pole masses (dashed) and a quenched 
(yellow) and unquenched (purple) LQCD calculations. The shaded bands reflect 
the theoretical uncertainties and are shown in the q 2 ranges for which calculations 
are reported. ([7]). 



Table 6: Representative predictions for flavor changing neutral current charm 
decays and experimental upper limits. 



Decay Mode 


LD (xl0 b ) 


MSSMR (xl0 b ) 


Experiment (90 % CL ul xl0 b ) 


D + -> vr+e+e- 


2.0 ® 


0.21(63) - 2.0(p|) 


7.4® 








11.2 dSZD 


D + -> ir + fj, + fj,~ 


1.9 


6.5(jMD - 15ff62D 


24.4([67]) 








3.9(168]) 


D+ -> vr + /i+e + 





30(62]) 


10.8(67) 


D u -> e+e" 


1.0 x 10 -1Y 


1.0 x 10" 4 


1.2(|69]) 




3.0 x 10~ 7 


3.5 


1.3© 








2.5 (70) 




1.8 


5.1 


ioo.o (HQ 
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Figure 4: The dilepton mass distribution for D + — > ir + e + e (normalized to Td+) 
in the MSSM with nonuniversal soft breaking effects. The solid line is the SM. (I) 
M- g = Mg = 250 GeV; (II) M g = 2M~ q = 500 GeV; (III) M g = M~ q = 1000 GeV; 
(IV) M g = (1/2) Mg = 250 GeV. Curves are from Ref. (El 




Figure 5: (a) The phase space corrected S-wave amplitude \S\ in the K + K~ 
system, shown as black points, from D° — > K + K~ir° decays. Red points are 
from an earlier analysis of D° — > K~ K + K$ decays. The solid blue curve is the 
line shape for /o(980) and the dotted curve for ao(980). In (b) are black points 
for |P| from D° — > K + K~it° decays. The solid blue curve is the </> line shape. 
The figure is taken from Ref. (I119p . 
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Table 7: Recent absolute branching fraction data. For CLEO-c results, the un- 
certainty due to radiative corrections has been absorbed into the systematic un- 
certainty. 



Mode 




Absolute BF (%) 


PDG(33| (%) 


D° —* 




4.007 ± 0.037 ± 0.070 ([86]) 


3.82 ±0.07 






3.891 ± 0.035 ± 0.069 ((HSJ) 




D° -> 


if-7T+7r 


14.57 ±01.2 ±0.069 ([82]) 




D° — > 


if - 7T + 7r + 7T — 


8.30 ±0.07 ±0.38 (82) 




D+ - 


► K+7T+7T+ 


9.14 ±0.10 ±0.17 ([82D 


9.51 ± 0.34 


D+ - 


► X+TT+TT+TT 


5.98 ±0.08 ±0.01 ([82D 


6.00 ± 0.28 


D+ - 




1.539 ± 0.022 ± 0.038 (El 


1.47 ±0.06 


D+ - 




7.05 ± 0.09 ± 0.25 ([82]) 


7.0 ±0.5 


D+ - 




3.149 ±0.046 ±0.096 ([82]) 


3.11 ±0.21 


D+ - 


► #+#-71-+ 


0.935 ± 0.017 ± 0.024 ([82]) 


1.0 ± 0.04 


D s -> 




5.50 ±0.23 ±0.16 ([83]) 


5.3 ±0.8 






4.0 ±0.4 ±0.4 flgZD 




D s — > 


K S K + 


1.49 ±0.07 ±0.05 dEHD 


2.2 ±0.45 


D s — > 


K+K-tt+it 


5.62 ±0.33 ±0.51 ([83]) 








1.11 ±0.07 ±0.04 ([M]) 


1.22 ±0.23 


D s -► 




1.47 ±0.12 ±0.14 ([83]) 


2.11 ±0.35 




7T + 7/ / 


4.02 ±0.27 ±0.30 (83) 


4.7 ±0.7 




Figure 6: (a) Decay times for D° — > K + tt~ decays. The range includes negative 
values arising from measurement uncertainties. The solid curve is the result of 
the fit described in the text allowing for mixing and the dashed curve assumes 
no mixing (x' 2 = y' = 0). Contributions from various types of background are 
indicated, (b) Residuals of data bins from the fit with no mixing. The solid 
curve is the fit with mixing, (c) Ratio of WS to RS decays in time slices. The 
approximately linear form expected from Eq. [42] is evident in the data. Figures 
are from Ref. (1124p (BaBar collaboration). 
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Table 8: Constraints on NP 


models from D° mixing. 


Model 


A " 4- /"I 4- ' 4- 

Approximate Constraint 


fourth Generation 


\Vub'Vcb>\ - my < 0.5 (GeV) 


<y = —l/o bmglet Ujuark 


s 2 ■ m s < 0.27 (GeV) 


C,/ = +2/6 binglet w;uark 


1 \ 1 ^- o A i n — 4 

\Kc\ < 2.4 • 10 


Little Higgs 


Iree: bee entry lor Q = —1/6 bmglet Quark 




Box: Region of parameter space can reach obsei 


Generic Z 


Mz'/C > 2.2 • 10 TeV 


T? *1 O 4- " 

family symmetries 


mi/ j > 1.2 ■ 10 lev (with m\/m2 = 0.5) 


Left-Right bymmetric 


No constraint 


Alternate Lett-Right symmetric 


Mr > 1.2 TeV (m^ = 0.5 TeV) 




/A / \ / 7\ T ~^ r\ A rxn X 7" — 1 

(Am/rriDj/Mji > 0.4 TeV 1 


Vector Leptoquark Bosons 


MyLQ > 55(App/0.l) TeV 


flavor Conserving 1 wo-Higgs-Doublet 


No constraint 


Flavor Changing Neutral Higgs 


m H /C > 2.4 • 10 3 TeV 


1C JNeutral Higgs (Cheng-bher ansatzj 


/ 1 A I ^ rr\r\ /"I ~\ j 

mul \A UC \ > b00 GeV 


Scalar Leptoquark Bosons 


See entry for RPV SUSY 


XT" 1 

Higgsiess 


TI IT -1 /"A _~\ rx^ X 7" 

M > 100 TeV 


X T 1 xn j_ T~\ " 

Universal Lxtra Dimensions 


JNo constraint 


O 1 " j_ X~i 

bplit lermion 


M/|Ay > (6 • I0 z GeV) 


Warped Geometries 


ivii > o.o lev 


IVTiniiTial Si rnersvin metric Standard 


I f^nh rriI < 3 5 • 10~ 2 for m ~ l TeV 




I0%)ll,hii| < -25 for m ~ l TeV 


Supersymmetric Alignment 


m > 2 TeV 


Supersymmetry with RPV 


A'mA'nfc/m^^l.S-lO^/lOO GeV 


Split Supersymmetry 


No constraint 



Table 9: Mixing and CPV parameters from D — > K + tt decays. 
Fit Type Parameter Fit Results / 10 -3 



BABAR(l24) CDF^Ml) Belle(l49) 



No CPV or mixing 




3.53 


± 


0.09 








3.77 ± 


0.01 




R D 


3.03 


± 


0.19 


3.04 


± 


0.55 


3.64 ± 


0.17 


No CPV 


x' 2 


-0.22 


± 


0.37 


-0.12 


± 


0.35 


0.18±°;| 




y' 


9.7 


± 


5.4 


8.5 


± 


7.6 


0-6^ 






R D 


3.03 


± 


0.19 














« D 


-21 


± 


54 








23 ± 


47 




a M 














670 ± 


1200 


CPV allowed 


x a+ 


-0.24 


± 


0.52 








< 0.72 






x> 2 ~ 


-0.20 


± 


0.50 














y' + 


9.8 


± 


7.8 








-28 < y' < 21 




y'- 


9.6 


± 


7.5 
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Figure 7: Dalitz plot of the decay D — > Kstt + tt showing squared invariant mass 
of i^57r + vs. that of Kstt~ for data from Ref. (|75p (Belle collaboration). 



Table 10: Mixing and CP V parameters from D — > /fs7r + 7r _ decays. The first 
uncertainty is statistical and the second systematic. The third is due to uncertain- 
ties in the isobar structure assumed in the model for the Dalitz plot distribution. 



Fit Type 


Parameter 


Fit Result ! 


35% C.L. Interval 


No CP 


x(%) 


0.80 ±0.29 t'm -14 


(0.0,1-6) 


Violation 


y(%) 


0.33 ± 0.24 +: °* 


(-0.34,0.96) 




x{%) 


0.81 ±0.30 ±f 7 tf 6 


|a?| < 1.6 


CPViol. 


y(%) 


0.37 ±0.25 -m 


\y\ < 1.04 


Allowed 


\Q/P\ 
arg q/p 


0.86+-^ + -° l ± 0.08 

r_u +16+5 +2 )° 

{ i4 -18 -3 -4) 




No direct 


\q/p\ 


0.951;^ 




CP viol. 


arg q/p 


C-2±S)° 
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Figure 8: (a) Likelihood contours projected onto (x D ,y D ) plane from 8 param- 
eter fit to 26 mixing observables for which data exists, (b) Projection onto the 
(\p/q\,axgp/q) plane for this fit. 



Table 11: Asymmetries in decays of D mesons to various final states /. The 
asymmetry is defined in the text. 



/ Acp{%) f A OP (%) 

D° decays: D + decays: 



Cabibbo Favored 






K-TT+ 


-0.4 ±1.0 




-0.9 + 0.9 


K s ir° 


+0.1 ± 1.3 


K S TT + 1T 


+0.3 + 0.9 




+0.2 ±0.9 




+0.1 ± 1.3 




+0.7 + 1.0 


K~TT + TT + 
K-TT + TT + TT° 


-0.5 + 1.0 
+1.0 + 1.3 


Cabibbo Suppressed 






K-K+ 


+0.15 + 0.34 


K S K+ 


+7.1 + 6.2 


K S K S 


-2.3 + 1.9 


K+K-TT+ 


+0.6 + 0.8 




+0.02 + 0.51 


K+K-TT+ 


+0.6 + 0.8 


TT°TT° 


+0.1 + 4.8 




-1.7 + 4.2 


7T~7r + 7r 


+0.1 + 5 


KsK + ir~ir + 


-4.2 + 6.8 


K-K + it--k + 


-8.2 + 7.3 






i£s7r + 7r~ 


-0-911? 






Doubly Cabibbo Suppressed 


Cabibbo Suppressed D s 


K+-K- 


-0.8 + 3.1 


D s -» K+r] 


-20 + 18 


K+TT-TT 


-0.1 + 5.2 


D s -> K+rj' 


-17 + 37 




-1.8 + 4.4 


D s -> K S TT+ 


27+11 



D s -> K + tt° 2 + 29 



